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ABSTRACT 


An  experimental  measurement  program  has  been  undertaken  to 
determine  the  parameters  for  propagation  in  the  40  and  70  Hz  range.  A 
transmitter  in  Wisconsin  (WTF)  transmitted  simple  sinusoidal  signals  for 
two  eight-hour  periods  per  day  for  20  days.  The  radiated  power  was  1/4  watt 
at  45  Hz  and  1/2  watt  at  75  Hz.  Receiver  sites  were  located  in  Utah,  Nova 
Scotia  and  Hawaii. 

Subject  to  the  constraints  detailed  in  the  text,  estimates  of  the 
attenuation  rates  and  the  mode  excitation  factors  have  been  determined  with 
high  accuracy.  The  daytime  attenuation  rates  were  found  to  be  higher  than 
those  estimated  previously  on  the  basis  of  theoretical  calculations.  The 
nighttime  excitation  factor  was  also  found  to  be  a  few  dB  smaller  than  pre¬ 
viously  expected.  A  careful  analysis  shows  that  directional  differences  in 
the  attenuation  rates  (  I  )  are  1 ess  than  0.  2  dB  |Mm  for  both 

the  40  and  70  Hz  range. 


Accepted  for  the  Air  Force 
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NOTATION 
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true  attenuation  rate  (dB/Mm)  of  only  propagating  mode 
true  attenuation  rate  (nepers/Mm) 

approximation  for  a  when  total  field  strengths  are  used  rather 
than  the  direct  field  components 

calibration  frequency  (Hz) 

signal  frequency  (Hz) 

excitation  factor  of  only  propagating  mode  (£  =  1  /YiJa^(c/v  ^)  ) 
phase  velocity  of  only  propagating  mode 

atmospheric  noise  (dBH/^Hz)  referred  to  an  equivalent  H  field 

<P 

voltage  gain  in  system 
calibration  monitor  voltage  (mv.  ) 

voltage  associated  with  the  received  field  strength  measured 
at  the  input  to  correlators. 

voltage  associated  with  the  injected  calibration  signal  measured 
at  the  input  to  correlators. 

voltage  at  the  antenna  output  associated  with  the  received  H 

9 

injected  calibration  signal  voltage  measured  across  the  input 
transformer 

voltage  attenuation  ratio  (associated  with  noise  processing 
circuit) 

distance  from  transmitter  to  receiver 
earth's  radius  -  6.378  Mm 
magnitude  of  a  quantity 
expectation 
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SUBSCRIPTS 


c 

cosine  component 

s 

sine  component 

c 

calibration  frequency 

s 

signal  frequency 

H 

Hawaii 

U 

Utah 

NS 

Nova  Scotia 

u 

unbiased  estimator 

b 

biased  estimator 

& 

direct  component 

fi 

round-the-world  component 

T 

total  field  component  (or  associated  with  gain  of  cascaded  tape 
recorders ) 

I 

integrated  (usually  refers  to  the  integrated  noise) 

i»  j 

indices 

D 

daytime 

N 

nighttime 

C.  M. 

calibration  monitor 

EW  (WE) 

propagation  direction  -  from  East  to  West  (from  West  to  East) 

m 

mid  or  center  frequency  in  a  selected  group  of  frequencies 

A 

(A2) 

SUPERSCRIPTS 

2 

estimator  for  the  true  value  of  A 

(A2) 

single  average 

I2 

double  average 
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I. 


INTRODUCTION 


For  the  period  15  March  through  9  April  1971,  the  test  facility- 
transmitter  at  Clam  Lake,  Wisconsin  was  scheduled  to  transmit  a  simple 
sinusoidal  signal  for  two  eight  hour  periods  per  day.  The  EW  antenna  con¬ 
sists  of  a  horizontal  dipole,  center  fed  with  the  ends  grounded.  The  dipole 
is  a  single  wire  transmission  line  hung  on  conventional  telephone  poles.  The 
total  length  of  the  EW  antenna  is  22.  5  km  and  the  line  current  was  300  A  (rms). 
The  electrical  axis  of  the  EW  antenna  is  approximately  116°  S.  E.  (There  is 
in  addition  a  similar  N.  S.  antenna  which  was  not  utilized  in  these  tests.  ) 

The  transmissions  were  specified  to  be  the  odd  frequencies  from  41  to  49  Hz 
and  from  73  to  77  Hz.  Receiver  sites  were  located  at  Wendover  -  Utah, 

Maitland  Bridge-Nova  Scotia  and  Pohakuloa  Test  Range-Hawaii  (Fig.  1). 

The  measurements  were  completed  without  significant  equipment  failures  or 
data  loss  for  the  duration  of  the  tests.  About  120  analog  magnetic  tapes  (40/ 
site),  each  with  eight  hours  of  recorded  signal  were  collected.  These 
tapes  have  been  processed  through  the  laboratory  playback  facility  and  the 
results  are  compiled  in  this  report. 

It  must  be  realized  at  the  outset  that  the  transmitted  signal  levels 
were  extremely  low.  For  example,  the  total  radiated  power  at  45  Hz 
was  only  l/4  watt  and  at  75  Hz  only  l/2  watt.  Since  the  Hawaii  site  is 
6.  53  Mm  away, to  obtain  good  estimates  of  the  signal  strength  one  needs  to 
perform  the  detection  with  equivalent  bandwidths  of  the  order  of  a  few  ten 
thousandths  of  a  Hz  (i.  e.  ,  integration  times  of  the  order  of  hours). 
Furthermore,  to  calculate  good  estimates  of  the  attenuation 
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Fig.  1.  Receiver  site  locations  (past,  present  and  future). 
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rate  a  it  is  necessary  to  average  several  such  (properly  weighted)  estimates 
at  each  receiver  site. 


The  propagation  parameters  of  interest  are  the  attenuation  rate  a  and 
a  term  denoted  as  the  excitation  factor,  1  /(h  *Joe  c  /v  (See  Eq.  1,  Section 

II  for  the  definition  of  terms).  These  two  quantities  vary  with  frequency, 
time  and  possibly  with  propagation  direction.  It  is  an  a  priori  hypothesis 
based  on  our  previous  limited  experimentation  that  both  of  these  quantities 
differ  significantly  from  all  day  to  all  night  path  conditions  but  remain  fairly 
constant  during  either  period.  The  propagation  direction  dependence  for 
30  <  f  <  80  Hz  is  believed  to  be  so  slight  that  an  extremely  precise  experi¬ 
ment  would  be  necessary  to  resolve  directional  differences. 

A  number  of  comprehensive  reviews  on  the  status  of  measurements 
and  theoretical  estimates  of  the  propagation  parameters  at  the  lower  ELF 
range  have  been  written  [  1,2,3]  .  Basically,  the  shortcomings  of  past 
parameter  estimates  fall  into  two  classes: 

(1)  Theoretical  estimates  depend  on  an  accurate  representation  of 
the  conductivity  profiles  of  the  D,  E  and  lower  F  regions  (including  effects  of 
ions  and  the  earth's  magnetic  field).  Small  differences  in  the  chosen  models 
lead  to  estimates  of  the  propagation  parameters  so  varied  that  an  ELF 
communications  system  could  not  be  sized  cn  theoretical  predictions 
alone . 


(2)  The  experimental  measurements  are  based  almost  entirely  on 
spectral  decomposition  of  individual  large  lightening  pulses.  Using  two 
receiver  sites  and  an  equation  similar  to  equation  1  (but  modified  to  account 
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for  a  vertical  rather  than  a  horizontal  dipole  source)  the  ratio  of  field  strengths 
at  the  two  sites  is  taken  in  order  to  estimate  o'(f).  This  method  is  rather 
unsuccessful  for  the  following  reasons.  First  the  point  of  origin  of  the 
particular  lightening  stroke  analyzed  is  imperfectly  known.  In  addition  it  is 
not  known  whether  the  stroke  orientation  is  cloud  to  cloud  (horizontal)  or 
cloud  to  ground  (vertical).  The  geometrical  configuration  of  the  receivers 
with  respect  to  the  lightening  stroke  may  be  such  that  small  uncertainties  in 
either  the  stroke's  position  or  in  the  signal  strength  estimate  may  give  large 
differences  in  values  of  a .  Second,  for  information  concerning  propagation 

at  frequencies  less  than  100  Hertz,  one  needs  to  analyze  the  long 
tail  of  the  stroke  in  question  which  tends  to  be  corrupted  with  return  strokes 
from  the  same  burst  and  the  bursts  of  other  strokes.  Therefore,  the  S/N 
ratios  tend  to  be  extremely  unfavorable  for  f  <  100  Hz.  Third,  on  the  tail  of 
the  sferic  there  may  be  included  a  component  of  signal  which  has  traveled 
around  the  world  and  adds  to  the  direct  path  component.  This  will  give  an 
additional  source  of  error  in  estimation  of  a. 

It  should  be  noted  that  recent  sferic  measurements  (Hughes  and 
Thiessen[4]  )  have  suggested  that  the  diurnal  difference  in  attenuation  rates 
is  about  1.  5  dB/Mm  at  50  Hz  (daytime  attenuation  being  the  higher).  If  one 
assumes  a  nighttime  attenuation  rate  of  .  7  dB/Mm  then  a ^  would  be  of  the 
order  of  Z  .  Z  dB/Mm.  (An  earlier  paper  by  Hughes  [5]  suggested  that 
significant  differences  in  and  exist.  )  In  order  to  determine 

whether  the  first  results  of  Hughes  (i.e.  ,  ~  Z.Z  dB/Mm)  were  valid  we 

deployed  receivers  in  Hawaii  and  California  in  an  attempt  to  measure 
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[6]  ).  This  experiment  was  completed  in  July  1970.  The  results  from  the 
test  were  somewhat  ambiguous  since  only  the  California  data  could  be 
utilized,  making  this  a  one  station  measurement.  The  Hawaii  data  was  not 
useable  since  the  summer  noise  was  considerably  higher  than  expected. 
(The  transmitted  signal  strength  was  also  some  6  dB  lower  than  at  present.  ) 
As  a  rule  of  thumb  for  the  accuracies  that  were  needed  in  that  experiment  to 
measure  o',  a  S/N  ratio  of  at  least  15  dB  is  required.  Because  of  low 
transmitted  signal  and  high  noise  conditions, even  after  a  3.2  hour  coherent 
integration  the  S/N^  ratio  in  Hawaii  was  only  8  dB.  Consequently,  we  were 
forced  to  base  our  conclusions  on  a  one  station  (California)  measurement  and 
make  assumptions  concerning  the  mode  excitation  factor 


( 1  /(h  J  <j^  c /Vp^) )  f°r  day  and  night  conditions  in  order  to  deduce  a  probable 
range  of  values  for  the  attenuation  rate  for  day  and  night.  The  tentative 
result  based  on  the  hypothetical  excitation  factors  was  that  was  reasonably 
close  to  .75  dB/Mm  at  45  Hz  while  the  value  of  was  probably  closer  to 
2  dB/Mm.  This  high  value  of  was  not  expected  and  was  a  strong  factor 
in  planning  the  recent  propagation  tests.  The  recent  tests  (with  the  two 
station  -  Utah  and  Hawaii  sites)  have  led  to  a  different  interpretation  of  these 
California  results.  Specifically  one  of  the  main  results  of  this  experiment 
was  that  the  nighttime  mode  excitation  factor  is  about  2.3  dB  lower  than 
originally  predicted  and  that  the  nighttime  attenuation  rate  is  about  .9  dB/Mm 
(at  45  Hz).  Hence  the  low  nighttime  signal  strengths  in  California  were  due  to 
a  worse  excitation  factor  than  originally  estimated  rather  than  a  high  attenuation 
rate;  in  this  perspective  the  earlier  tests  and  the  present  tests  are  totally 
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compatible. 

This  report  has  a  dual  purpose.  First  it  is  a  compilation  of  the  re¬ 
sults  of  the  recent  experiment.  Second  it  is  intended  for  use  as  a  handbook 
in  that  it  describes  in  considerable  detail  the  receiver  design  considerations, 
the  mechanics  of  the  system  calibration  and  the  methods  by  which  the  data  is 
analyzed  and  estimates  of  a  and  the  excitation  factor  are  determined.  Since 
this  type  of  experiment  is  to  be  repeated  a  few  more  times  by  both  Lincoln 
Laboratory  and  Navy  laboratories  this  report  should  serve  as  a  useful 
reference  with  regard  to  standardization  of  experimental  technique. 

For  those  interested  only  in  the  detailed  results  of  the  recent  experi¬ 
mental  phase  it  is  probably  sufficient  to  read  the  main  text  (Sections  I 
through  VI).  The  conclusions  are  presented  in  Section  VI.  Appendices  A 
through  I  serve  as  reference  material. 

The  primary  goal  of  the  recent  experimental  phase  was  to  measure 
^EW  (a^enuari°n  rate  in  an  east  to  west  direction)  using  data  from  Hawaii 
and  Utah.  We  believe  we  have  determined  this  value  satisfactorily. 

A  second  goal  was  to  use  the  data  from  Utah  and  Nova  Scotia  to  look  for 

small  differences  (i.  e.  ,  ±  .  2  dB/Mm)  in  (  a-r “  0t,,rTP)-  We  were  able 

ili  W  W  ill 

to  achieve  this  goal  as  well. 

All  field  strengths  are  expressed  in  terms  of  H  (cp  denotes  a 

horizontal  component  normal  to  the  propagation  direction).  The  magnitude 

of  the  vertical  E  field  is  related  to  H  through 

cp 
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where  ^ 


=  377fi 


and  is  the  phase  velocity. 


II. 


THEORY 


The  equations  determining  the  ELF  electromagnetic  field  com¬ 
ponents  produced  by  a  specified  source  at  the  surface  of  a  spherical  earth 
and  confined  by  the  presence  of  a  concentric  and  vertically  inhomogeneous 
ionosphere  have  been  extensively  investigated  over  the  last  decade.  One 
basic  formulation  is  in  terms  of  a  zonal  harmonic  series.  The  use  of  this 
method  is  illustrated  in  a  series  of  calculations  performed  by  Johler  and 
Lewis  [7]  for  a  specified  ionospheric  model  and  frequency  of  operation.  The 
disadvantage  of  this  formulation  is  that  the  equations  (which  have  a  few 
hundred  significant  terms  in  the  series)  are  not  readily  amenable  to  confirma¬ 
tion  within  a  reasonable  experimental  program.  In  addition, because  of  the 
complexity  of  the  series  it  is  difficult  to  develop  a  physical  insight  into  the 
nature  of  the  propagation  and  hence  it  is  difficult  to  extrapolate  the  results  to 
other  frequencies  or  models  where  the  ionosphere  has  markedly  different 
characteristics  over  different  portions  of  the  propagation  path  (i.e.  ,  daytime 
versus  nighttime  effects).  A  more  useful  (and  physically  revealing) 
formulation  results  from  use  of  the  Watson  transformation  (Wait  [8,9]  )  which 
converts  the  series  of  zonal  harmonics  to  a  more  rapidly  converging  infinite 
series  of  radial  harmonics.  When  the  receiver  is  located  in  the  far  field 


2l rd 


>  1)  but  sufficiently  removed  from  the  region  of  the  transmitter's 


antipode,  and  the  attenuation  rate  is  greater  than  a  few  tenths  of  a  dB/Mm, 
the  series  of  radial  harmonics  can  be  replaced  to  first  order  by  the 
simple  expression  (for  a  horizontal  electric  dipole  transmitter  source 
at  the  earth's  surface)  representing  the  direct  ($)  contribution  to  the 
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H  field 
O 


where 

ILi 


ILf 

2-V 


cos  $  f- 


W<7, 


(—) 

eVph 


(1) 


j(k 


o 


X  e 


d-wfc  +  ff/4) 

Vph 


current  moment  of  the  transmitter 


$  angle  between  axis  of  horizontal  antenna  and  propagation  path 
direction 


' e 
h 


v 


ph 


V 


ground  conductivity  at  transmitter 
effective  height  of  ionosphere 
phase  velocity  for  the  only  propagating  mode 
120tt 


r  earth's  radius 

e 

kQ  2n/\ 

d  shorter  great  circle  distance  between  receiver  and  transmitter 

d  C  7 7  r 

e 

a?'  attenuation  rate  (nepers/m)  for  the  only  propagating  mode. 

For  greater  accuracy  one  must  add  to  the  field  shown  in  equation  1  another 


magnetic  field  component  (ft  component)  which  travels  around  the  world  along 


the  longer  segment  of  the  same  great  circle  path.  This  component  is  given  by 


(H  )B  =  Fe 

<p  rC 


-a'(2nre-d)  j[k0(^— )  (27rre-d)] 


(2) 


where 
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ILf  /Z75T  „  r 

=  ~TZ  V— 7 cos  $  1“ 


2r7 


X  e 
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}  r- 


*e  J  e 


ph 


—  r.  d- 

r  J  sm  — 


(3) 


■jM  -f). 


These  direct  and  round  the  world  components  add  at  each  point  to  give  a 
standing  wave  pattern  as  illustrated  in  Fig.  2.  It  can  be  shown  (Jones  [  10]  ) 
that  the  minimums  of  the  standing  wave  H  pattern  are  located  at  a  distance 
p'  from  the  antipode  where 

,  _  (4n+ 3 )  X 

p  ' 

and  the  standing  wave  length  is 


n  =  0,  1 ,  2,  ... 


(4) 


s  w 


—  I  j  ^  A 

^n+l  ^ n  2 (c/Vph)  ’ 


(5) 


It  should  be  made  clear  that  the  total  field  (H  )rn  as  expressed  by  the 

o  1 

sum  of  equations  1  and  2  was  derived  under  the  set  of  assumptions 


a)  The  ionosphere  is  variable  only  in  the  vertical  direction. 
No  lateral  inhomogeneities  (i.e.  ,  day-night  transitions)  are  allowed. 

b)  The  ionosphere  is  isotropic,  that  is,  the  earth's  magnetic 
field  is  assumed  to  be  zero.  This  is  a  reasonable  approximation  provided 
that  a  significant  amount  of  the  electromagnetic  energy  does  not  get  into  the 
altitude  region  where  the  mean  free  path  between  electron-neutral  collisions 
becomes  comparable  to  the  electron  gyro  radius.  During  daytime  conditions 
this  appears  to  be  the  case  but  for  nighttime  conditions  the  magnetic  field 
probably  has  some  minor  effect  on  the  directional  behavior  of  the  attenuation 
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Fig.  Z.  Standing  wave  pattern  in  an  idealized 
uniform  earth-ionosphere  cavity. 
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rate. 


c)  The  three  parameters  h,  a  and  v  ^  are  expected  to 
reflect  the  behavior  of  the  lower  ionosphere  (and  to  lesser  extent  the  variable 
ground  conductivity)  on  the  ELF  propagation.  Only  an  extensive  experimental 
program  will  show  if  these  parameters  are  sufficient  to  characterize  the  bulk 
of  the  propagation  data. 

d) 

The  form  of  equations  1  and  2  suggests  that  a  good  approximation  for 


<y,  h  and  v  ^  are  all  frequency  dependent. 


the  total  field  strength  (H  in  a  more  realistic  model  for  the  earth- 

<P  1 

ionosphere  cavity  where  the  changes  in  a  and  v  ^  can  be  both  locally  and 
directionally  varying  is  (Fig,  3) 

-J r  "yr)dr  Jko  /  (c/vphVdr 

(H  )  =  F'e  °  e  ° 

<p  a 


and 


2  77  r  -d 
e 


277  r  -d 
e 


P  oi'  (r)dr  jk  f  (c/v 

Jo  9  J  oJo  ' 


(H  )  =  +  F  'e 

<p  ft 


ph  9 


dr 


(6) 


(7) 


and 


(H  )  =  (H  )  +  (H  ) 

05  (0  (/5  ft 


(8) 


F'  differs  from  F  in  the  following  way:  h*/fc:/v  is  replaced  by  an  effective 

number  (h^t/v  ^)ce)  which  can  only  be  determined  from  experimental  data. 
The  primary  objective  of  this  measurement  program  was  to  estimate  the 

magnitude  of  cv  (r).  With  CW  signals  of  the  type  radiated  by  the  WTF,  one  can 

& 

by  experiment  estimate  the  magnitude  of  the  total  horizontal  magnetic  field 
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Fig,  3.  Pictorial  representation  of  the  standing  wave  pattern 
in  a  realistic  model  of  the  earth-ionosphere  cavity. 
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(H  )_.  The  difficulty  comes  in  relating  the  estimates  of  (H  )  and  the 
CP  1  (D  1 

function  a  (r),  The  first  hypothesis  then  is  that  the  exponential  terms  in 

$ 

equations  6  and  7  are  approximately  constant  over  the  estimate's  determina¬ 
tion  time.  Accordingly  since  the  estimation  period  is  of  the  order  of  a  few 
hours  and  the  direct  path  length  is  at  least  a  few  megameters  then 


dr>  d^ 

'  «  (r.tldr-Cj,  J  (c/v  ,)  d 


r  1 


r  1 


i  /  -r  -  C, 

ph  &  3 


(9) 


and 


2  77  r  -d  , 
e  r  1 


2 77 r  -d 


J 


277 r  "d 
e  rZ 


aq  (r ,  t)dr  C2  ,  J 


e  r  1 


277  r  -d 
e  r  Z 


(c/v  ,  )  dr  ^  C 


for  all  t^  £  t  £  t^  where  and  t^  are  the  beginning  and  end  of  the  estimation 
period,  d^  and  d^^  denote  the  locations  of  two  receivers  on  the  same  great 
circle  path  as  the  transmitter.  From  prior  experiments  [6]  it  has  been 
determined  (as  far  as  poor  S/N  ratios  will  allow)  that  as  long  as  the  entire 
direct  path  segment  is  in  daytime  (or  nighttime)  conditions  the  received  signal 
strength  is  approximately  a  constant  value  over  a  period  of  several 
hours.  The  only  observable  signal  changes  come  when  the  sunrise  (or  sunset) 
terminator  intercepts  the  direct  path  segment.  Hence  for  the  remainder  of 
this  report  we  assume  that 


& 


for 

(tss 

*,d(£) 

(tSR 

lSR 

fcss 


(10) 
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and 


(c/vph(f)) 

(c/vpht£)) 


#.N 


(ID 


where  tcr)  and  tcc  are  the  appropriately  chosen  ionospheric  sunrise  and 
oK  bo 

sunset  times  and  26  is  an  appropriate  transition  time  interval.  D  and  N 

refer  to  daytime  and  nighttime  conditions. 

It  must  be  emphasized  that  the  conditions  expressed  in  equations  10  and 

11  are  very  strong  hypotheses  that  at  present  we  are  forced  to  make.  The 

reason  is  that  for  estimation  purposes  one  needs  an  integrated  signal  to  noise 

ratio  (S/Nj)  of  at  least  +  20  dB  if  we  are  to  determine  a  good  estimate  of 

(H  )_,.  However  at  a  typical  receiver  site  the  S/N  ratio  (as  measured  in  a 
<P  1 

1  Hz  bandwidth)  is  usually  about  -2  5  dB.  Hence  one  must  integrate  for  a  few 

hours  to  achieve  a  sufficient  value  of  S/Nj  for  the  purpose  of  estimating 

(H  with  a  reasonable  set  of  error  bounds.  In  order  to  deduce  useful 
O  T 


estimates  of  a.  (f)  from  the  (H  estimates  one  needs  the  error  bounds  on 
A  6  f  T 

(H  )rn  that  would  result  from  integrating  over  tens  of  hours.  It  is,  of  course, 

< 0  T 

necessary  that  the  signal  amplitude  and  phase  are  nearly  constant  over  this 

long  integration  time.  To  approximate  this  long  integration  time  we  take  a 

few  equispaced  samples  per  day  (or  night),  collect  samples  over  a  few  days 

and  incoherently  average  the  results.  It  must  be  emphasized  that  this 

technique  presupposes  that  ^^(f)  and  (c/v  N(D)  are  essen^a^y 

constant  from  day  to  day  (particularly  for  the  same  GMT  time  intervals). 

The  bulk  of  our  experimental  data  supports  this  last  assumption  (again  within 
— 

See  schedule  for  method  of  selecting  tgg  and  tg^  (Appendix  H). 


15 


the  limitations  of  the  generally  low  S/N  ratios). 


We  have  noted  that  the  measurements  performed  in  this  experiment 

yield  |H  (d)  |  and  not  |H  (d)  I  .  (Only  the  magnitude  is  measured 
1  <p  1  (0  » 

and  not  its  relative  phase.  )  Normally  |H^(d)  |  ^  is  measured  at 

a  few  receiver  sites  all  on  the  same  great -circle  path  with  the  transmitter. 

Since  we  do  not  know  the  magnitude  of  (H^(d)^  or  its  phase  relative  to 

(H  (d))  one  cannot  directly  deduce  (H  (d))  from  (H  (d))^.  However,  for 

<p  $  <p  J9  <p  1 

all  the  receiver  sites  where  we  have  been  or  plan  on  operating  it  is  true  that 

|H  (d)  |  >  >  |H  (d)  |  .  Hence  for  the  range  (1  <  d  £  10  Mm)  the  total  field 
1  o  $  (p 

standing  wave  pattern  appears  as  a  sinusoidal  pattern  with  an  amplitude  which 

increases  with  d  superimposed  on  the  direct  field  (see  Fig.  2).  Then,  if 

the  relative  phase  between  (H  )  and  (H  )  at  a  particular  site  were  random 

tp  (p  ft 

from  measurement  to  measurement  it  seems  obvious  that  an  appropriate 

average  of  a  series  of  |(H  )_  I  measurements  would  give  a  value  very  close 

(0  1 

to  the  magnitude  |(H  )  I.  ‘  If  this  were  done  at  two  receiver  sites  on  the  same 

(p  $ 

great  circle  path  (see  Section  IV  for  details)  one  could  then  compute  a  value 

for  a  (f)  which  is  relatively  free  of  errors  due  to  the  standing  wave  pattern. 

$ 

The  problem  reduces  to  providing  relative  phase  shifts  between  (H  )  and 

cp  $ 

(H  )  for  a  large  enough  set  of  measurements.  A  simple  way  to  change  phase 

O  ft 

is  to  change  frequency  slightly  (one  must  assume  that  the  attenuation  rate 
does  not  change  significantly  over  this  same  small  frequency  shift).  By 


The  probability  distribution  function  for  the  relative  phase  between  (H  ) 

(p  $ 

and  (H  )  must  however  be  reasonably  uniform  from  0  to  2flr. 
tO  ft 
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selecting  a  number  of  slightly  different  frequencies  a  proper  set  of  phase 
shifts  can  be  obtained  and  hence  averaging  the  total  fields  will  give  a  result 
close  to  the  direct  field  alone.  (The  theoretical  justification  for  this  pro¬ 
cedure  is  detailed  in  Appendix  G.  ) 
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III.  DESCRIPTION  OF  EXPERIMENTAL  APPARATUS  AND  METHOD  OF 
ANALYSIS 

Figures  4  and  5  are  simplified  block  diagrams  showing  the  field 
site  recording  equipment  and  the  laboratory  playback  facility.  A  detailed 
description  of  the  field  site  recording  equipment  is  presented  in  Appendix  A. 
This  section  provides  a  brief  discussion  of  how  the  receiver  processes 
the  received  signal  and  permits  a  simple  calibration  procedure. 

An  injected  calibration  signal  (whose  frequency  is  1  Hz  removed  from 
the  transmitted  frequency)  is  used  to  provide  a  continuous  gain  calibration 
of  the  receiver.  A  frequency  synthesizer  provides  the  calibration  signal  as 
well  as  the  phase  reference  signals  which  are  used  in  the  quadrature 
correlation  receivers.  The  calibration  procedure  is  discussed  more  fully 
in  Appendices  E  and  A.  Note  from  Fig.  4  that  both  a  narrowband  and  a 
wideband  channel  are  provided.  All  of  the  processing  reported  here  is  on  the 
narrowband  channel  data.  All  the  data  were  recorded  on  magnetic  tape 
and  processed  at  the  laboratory  playback  facility. 

The  magnetic  tapes  are  played  back  (see  Fig.  5)  with  a  speedup  factor 
of  32  so  as  to  reduce  processing  time  from  about  12  hours  per  tape  to  about 
23  minutes.  Each  of  the  channels  consists  of  a  pair  of  correlators  in  quadra¬ 
ture  followed  by  squarers  and  a  summer.  The  outputs  are  recorded  on  a 
chart  recorder.  (A  Hewlett-Packard  2  pen  X-Y  recorder  provides  good 
resolution.  )  The  units  marked  PAR  (Princeton  Applied  Research  Model  120) 
consist  basically  of  a  multiplier  followed  by  a  single  pole  low  pass  filter 
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Fig.  4.  Field  site  recording  equipment. 


r  ( t ,  n )  ( Input  Signal)  =  n  (U,t)  +  A  cos  (wst  +  <f>  +  8cos  (a>c  t  +  G)  ~Tl8 -6-13138  2 


o 
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+ 
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Fig.  5.  Playback  facility  block  diagram. 


(see  Fig.  E-l).  It  can  be  shown  that  the  equivalent  integration  time  is 
T  =  2(32)RC  where  RC  is  the  time  constant  of  the  low  pass  filter.  For  most 
of  our  processing  RC  =  80  sec  was  used  which  corresponds  to  T  —  1. 42  hours. 
Even  with  this  long  integration  time  the  signal  to  noise  ratio  at  Hawaii  was 
typically  only  10-20  dB.  This  low  value  of  S/N^  (low  for  purposes  of 
accurate  estimation  of  the  signal  level)  necessitates  a  very  careful  study  of 
the  noise.  As  noted  in  Appendix  A  the  noise  of  interest  is  distinctly  non- 
Gaussian  atmospheric  noise  which  completely  dominates  the  receiving  system 
noise.  The  input  signal  to  the  quadrature  correlator  receivers  can  be  represented 

r(t,  0)  -  n(t,  fi)  +  Acos  (2Trfgt  +  ^)  +  B  cos  (2Trf^t  +  A)  (12) 

where  (A,<n)  and  (B,  0)  are  the  amplitude  and  phase  angles  associated  with 

the  transmitted  and  calibration  signals  respectively.  n(t,  .  )  is  a  non-Gaussian 

noise  process  with  zero  mean.  Let  us  consider  the  transmitted  signal  channel 

only  since  the  results  for  the  calibration  signal  are  the  same  with  (B,  0) 

replacing  (A ,(p).  Let  the  outputs  of  the  PARs  at  time  T  be  denoted  by  the  random 

variables  L  and  L  .  It  is  an  easy  matter  to  show  that 
c  s  3 

?  T 

E[L^]  =  E[—  J  j:os go  t)  r  (t,  fi)dt]  =  A  cos  (p 
o 

2  T 

E[Lg]  =  E[^  J  fcin  a>gt)  r(t,  Q)dt]  =  A  sin  ^ 
o 

var  [LJ  =  var  [Lg]  =  Nq/T  =  Nj  (13) 

where  n(t,  .  )  was  assumed  to  be  a  white  noise  process  with  the  spectral  density 
given  by 


2  1 


-  00  <  f  <  +  00  t 


(14) 


N 

S  (f)  =  -f, 
n  Z 

Lincoln  Laboratory  wideband  ELF  noise  studies  [  11]  have  shown  the  noise 

to  be  nearly  'white'  over  the  bandwidths  used  in  this  set  of  measurements. 

N  is  termed  the  integrated  noise.  After  L  and  L  are  individually  squared 

sc 

and  then  summed,  the  result  is  another  random  variable  (denoted  by  A^) 
which  has  the  properties 

A 

E[A^]  =  a2  {1  +  ZNj/A2} 

and 


(15) 


var  [A^]  =  4(NI)2  {l  +  A2/^}  .  (16) 

2  2 

It  is  seen  that  the  estimator  A,  is  a  biased  estimator  of  A  with  the  bias 

b  - 

term  being  (2NJ.  After  the  determination  of  Ni  the  bias  term  can  be  sub- 
1  1  /s 

2  2 

tracted  off  to  give  the  unbiased  estimator  of  A  which  we  denote  by  A  to 

/N  /\  U 

2  2  2 

differentiate  it  from  the  biased  estimate  A^  (i.e.  ,  A^  =  A^  -  2Nj).  Note 
that  this  is  one  of  the  reasons  why  must  be  calculated.  See  Appendix  B 
for  an  explanation  on  how  is  calculated. 

We  will  now  outline  the  steps  one  must  take  in  order  to  obtain  a  single 
unbiased  estimate  of  input  field  strength. 

XN  /N 

2  2 

1)  At  a  specified  sample  time,  L,  values  of  A^(t^)  and  B  (L)  are 


read  off  the  calibrated  X-Y  chart  recorder  graph. 


2)  Using  these  two  values  plus  the  calibration  constants  measured 

✓ 

2 

in  the  field  the  biased  estimator  of  |H  L  is  computed  with  the  use  of 


< D 


equation  (e-8)  (Appendix  E). 


3)  The  equivalent  integrated  input  noise  (H^(L))j,  is  computed. 
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First  is  determined  as  shown  in  Appendix  B.  The  variance 

computed  from  the  Varian  620i  program  is  related  to  the  equivalent  input  H 

<P 

noise  (H^.  dBH/ ^Hz )  by  equation  b-19.  The  noise  typically  is  analyzed  for 

a  period  of  an  hour  centered  about  the  sampling  time,  t. ,  of  the  estimator 

/N  1 

A^U).  The  equivalent  integration  time  T  (=  2(32)RC)  is  determined  by  the 
quality  of  the  data  one  is  analyzing.  For  this  set  of  experiments  we  have 
chosen  RC  =  80  sec  which  sets  T  =  1.42  hours.  Then 


by 


4) 


(Hn}i  =  ((hn}/t  •  <17> 

The  unbiased  estimator  for  the  input  field  strength  is  then  given 


IH0(fS»u  =  [H„(fS»b-aHN>fS»/T- 


18) 


Figures  8  to  15  show  typical  X-Y  recordings  from  Utah  and  Hawaii.  The 
A2  ^2 

traces  are  and  B  versus  time  of  day.  The  samples  are  separated  far 

enough  in  time  that  they  are  independent  samples. 

The  question  of  'error  bars  *  on  the  sample  estimates  is  particularly 

relevant  since  it  is  important  to  be  able  to  differentiate  between  level  changes 

due  to  changes  in  signal  amplitude  or  phase  and  variations  due  to  the  normal 

noise  behavior.  This  is  particularly  important  since  we  have  made  the 

assumption  that  the  received  signal  level  is  approximately  constant  over  the 

recording  period  and  yet  we  commonly  see  changes  in  A^  over  this  period. 

A  b 

It  is  important  to  determine  whether  the  changes  in  A£  over  the  recording 
period  are  entirely  consistent  with  the  integrated  signal  to  noise  ratio  alone. 
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For  each  estimate  of  signal  strength  and  integrated  noise  one  can  compute 
a  sample  S/N  ratio 


<A2<Vfs»b 

No(,i'fS)/T 


b,T 


<ti’fS)/T 


(19) 


A 

The  quantities  (A‘"(T,fg))^  and  No(h,fg)/T  refer  to  the  output  of  the  receiver 
^2  Z 

whereas  (H  (t.,fc)),  and  H,T(t.,fc)/T  refer  to  equivalent  magnetic  field 
1  o  b9  1  lNib 

strengths  at  the  antenna  inputs.  Throughout  the  text  we  often  use  the  first 
pair  of  quantities  since  they  are  obtained  first  in  the  analysis.  This 
particular  sample  S/N^  ratio  then  uniquely  specifies  a  confidence  interval 
(for  a  specified  confidence  level  ^).  Since  the  uncertainty  in  is  negligible 
compared  to  the  uncertainties  in  S/N^  due  to  noise,  we  assume  that  is  in 
fact  a  known  quantity  as  far  as  determining  the  confidence  interval.  Appendix 
C  describes  in  detail  the  means  of  determining  the  confidence  interval 
(A^(e),  A^t(^))  for  each  field  strength  sample.  The  interpretation  is  that 
the  confidence  interval  (itself  a  random  interval)  includes  the  true  value  of 
A^  with  a  probability  (1  -  x  100$. 


The  integrated  signal  to  noise  ratio  obtained  for  a  single  sample  may 

be  satisfactory  for  an  estimate  of  |H  L.  However,  it  is  difficult  to  employ 

<p  1 

simple  H  estimates  from  two  stations  to  measure  an  effective 

<P 

attenuation  rate  Q'e££*  Let  us  assume  as  an  example  two  stations  separated 
by  4  Mm  with  the  stations  having  S/N  ratios  of  20  dB  and  30  dB.  (These 
numbers  are  greater  than  those  exhibited  in  Utah  and  Hawaii.  )  The  80$ 
confidence  limits  are  computed  to  be  (20  +  1.0  dB,  20  -  1.25  dB)  and 
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(30  ±  .  35  dB)  respectively.  Thus  the  noise  in  itself  leads  to  a  spread  in  &e££ 
estimates  of  about  ±  .  35dB/Mm.  The  obvious  step  then  is  to  incoherently 
average  successive  samples  where  approximately  the  same  ionospheric 
conditions  occur  and  if  this  is  not  enough  one  can  then  average  from  day  to 
day.  It  seems  obvious  that  one  should  also  weight  more  heavily  those  data 
where  the  noise  is  low.  One  appropriate  weighting  would  be  to  weight  each 
sample  inversely  by  its  measured  noise  sample.  At  the  frequency  T  the 
average  is  then 


M  (A.(f.)) 
v  v  r;u 


— ? —  i  nt  . 

<A2(f:))  = 


J 


i=lNI,i 


(20) 


or  the  equivalent  in  terms  of  H 


<P 


V  - _ o  I  u>  1?  1 


|H2  (f:) 


(D  j  'u,T  M 

L 


i=1(HN>I.i 


(21) 


i=1  (HN^I,i 


where  i  denotes  the  ith  sample  out  of  the  M  total  samples  taken  at  frequency 


fy  The  average  (A  (f.))  is  a  random  variable  and  Appendix  D  outlines  how 
to  specify  confidence  limits  for  this  random  variable. 

By  the  above  process,  at  each  frequency  f \ ,  a  weighted  average  value 
for  the  square  of  the  total  field  strength  is  calculated.  However,  what  is 
desired  is  not  the  total  field  (which  includes  the  round-the-world  component) 
but  rather  a  good  estimate  of  the  direct  component  of  the  field.  In  order 
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to  average  out  the  around- the -world  signal  the  weighted  total 


fields,  H^(f.)  are  in  turn  averaged  over  the  selected  frequencies  to  give  an 

o  J 

approximation  for  the  direct  field  contribution  alone  at  the  mid-frequency, 
f  .  The  second  averaging  process  result  in 


I  9  ,  2N+1  |H  (f.)  I  „ 

|H^(f  )  I  =  f^  - - -  — cfi — J — 2L».T  . 

'Vm^uJ  m(2N+l)  .  f2 

J  j 

2 

The  inverse  weighting  with  f  is  to  reduce  the  frequency  dependence  (see 

equation  1).  It  must  be  remembered  that  over  the  band  (f^ ,  .  .  .  ,  f^N+1^  >a  *  h 

and  c/v  ^  were  assumed  to  be  constant  with  frequency,  f^  is  the 

midpoint  frequency  in  this  band.  Also  explained  in  Appendix  D  is  the  means 

2 

of  calculating  confidence  limits  on  the  random  variable  |H  (f  )  I 

1  ip  m  lu,i 


(22) 
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IV. 


CALCULATION  OF  THE  ATTENUATION  RATE  AND  EXCITATION 
FACTOR 

Let  us  assume  that  we  have  computed  unbiased  estimates  for  the 


direct  field  strength  |H  I  at  two  receiver  sites  on  a  great  circle  which 
6  1  cp  jG>  u  6 

also  intercepts  the  transmitter  site.  Denote  the  receiver  site  locations  by 
x  and  y.  The  magnitude  of  the  field  strength  is  related  to  the  attenuation 
rate  by  equation  1  which  can  be  rewritten 

T77—  ■  ' 


H 


(f>  u 


ILf 


-Of 


cos  $  f- 


ph 


■)Mr? 


sm 


d T* 


}e 


•  (23) 


Let  us  consider  the  case  where  the  two  receivers  are  on  the  same  side  of  the 
transmitter.  The  above  equation  can  then  be  used  in  two  distinct  ways. 

First,  for  estimates  made  for  the  same  time  period  at  two  separated  sites 
we  have 


|H  I  /sin  d  /r  -a'  (d  -d  ) 

V &u'?-  =  f  J  .  -  -e  )  e  0  x  y  . 

i  sin  Vre 


Or 


ip  'jfi.u.y 


.  sin(d  /r  )  |H  | 

%  *  d^d~  tl01og  sin)///)  t201°8  =a/'U,X)  dB/Mm.  (24) 

y  X  y  '  lHAu,y 


It  is  evident  from  this  equation  that  the  excitation  factor  p  =  1  / (h^g^  (c /v  ) 

cancels  out  in  the  ratio  since  the  estimates  were  taken  for  the  same  time 

periods  for  both  sites.  In  this  manner  one  can  then  determine  a  ^  and 

_  _  D 

ot  .  f.T  by  using  the  pairs  of  estimates  f  |H  I  |H  I  _  }  and 

j£>.  N  7  6  ^  c  1  (p'fi.u.x.D  I  ^l^.u.y.D-' 

f  |H  I  ,  |H  I  }.  The  restrictions  on  a  are  those  referred  to 

1  <0  'jB,  u,  x,  N  1  1  u,  y ,  N  J  jB 
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in  the  text  associated  with  equation  10.  After  a  has  been  determined  the 

$ 

excitation  factor  (g,  =  1  / (h/s/o^T^Tv^T )  can  be  determined  from  equation  23 
since  all  the  other  terms  are  known. 


As  a  result  of  the  high  accuracy  to  which  |H  must  be  known  there 

may  be  a  residual  bias  error  in  calibration  which,  although  slight,  is 

significant  in  light  of  the  high  accuracies  we  require  in  the  |H  1^  measure- 

<p  1 

ments.  If  either  or  both  of  the  two  sites  has  a  consistent  bias  error  in  its 

calibrations,  this  error  can  show  up  as  a  potentially  large  error  in  the 

estimate  of  q^.  (i.e.  A  1  dB  error  in  estimating  a  field  strength  is  not 

usually  considered  serious  whereas  an  error  of  say  0.3  dB/Mm  in  the 

estimate  of  oi  can  have  profound  effects  in  sizing  a  global  communications 

system.  )  However,  we  can  obtain  good  relative  information  on  a  and  p 

J9 

even  when  a  bias  error  in  calibration  for  a  particular  site  is  significant.  We 
choose  the  signal  estimates  at  two  different  times  (a  day  estimate  and  a 
nighttime  estimate)  for  a  single  station.  The  following  ratio  is  then  computed. 


IH 


(p'£,u,x,  D  


fVc/v  ,  } 


ph  JN 


(a '  ta  -  a'  .  ,T)d 

jS,  D  jQ,  N'  x 


lH<pl,0,u,x,N  fVc/vph  }D 
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£,  D  N  x 


(25) 


Note  that  if  a  calibration  bias  exists  and  is  constant  from  day  to  night  then 
this  bias  will  cancel  in  the  above  equation.  This  can  be  rewritten  as 


H 


20  log 


jZl 


IH 


p 

iQ,  u,  x,  D  ^  ^  ,  D 

r -  =  20  l0g  7^  ‘  (“«,D  -  n fj,  N  dx* 


(26) 


(5  1  jS,  u,  x,  N 
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We  can  apply  the  same  sort  of  reasoning  to  the  receiver  at  location  and 
get 


2  0  log 


|H  I  „ 

lH  I  Kt 

1  ^  u,  y,  N 


20106  e^-S,D 


a 


fl.N 


It  is  evident  that  equations  26  and  27  are  two  linear  equations  for  the  two 

P  D 

unknowns  (a  ~  -  oi  ,T)  and  20  log  - .  The  solutions  are 

jfi.D  ft,  N  &  P.„ 


(27) 


(a  D  -  or  N>  =  H  1  .  f  20  log  ^alAUlryrD 
fl,D  *,N  dx-d  iH  | 

1  co  Jfi»  u*  y*  n 


2  0  log 


co  ft*  u»  y»  n 

.  20  log  Sru'x'D  J 

^<p  u,x,  N 

T2  =  d“7d-  [dx(2  0  1og  ^L^ytP) 

PN  x  y  |H  I  M 

1  cp'jB,  u.y.N 


|H  I  „ 

-d  (20  log  -=g  X>  — 

y  iu  i 


(28) 


(29) 


Thus  one  can  compute  the  difference  (cv  ~  -  a  -T)  and  the  excitation  ratio 

jR,  JJ  ft,  JN 

rd/rn  even  if  there  is  a  consistent  bias  uncertainty  at  either  or  both  of  the 
sites . 


It  should  be  clear  that  the  preceding  discussion  outlines  a  method  of 

determining  the  attenuation  rate  a  in  any  direction.  Since  in  this  series 

J0 

of  tests  two  of  our  receivers  are  in  Utah  and  Hawaii  we  can  determine 
estimates  of  aEW>D»  «EW)N>  &D  and  £N- 

We  should  also  like  to  investigate  any  anisotrophy  in  a  which  is  caused 
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by  the  earth's  magnetic  field  Bq  (i.e.  ,  the  difference  o -  ^WE^*  This 
can  be  investigated  by  spacing  two  receivers  equidistant  from  the  transmitter 
and  on  opposite  sides  of  the  same  great  circle  path  (d  =  d  =  d).  Moreover, 

x  y 

estimates  are  made  when  both  receivers  are  under  either  day  or  under  night 
conditions.  Since  the  receivers  are  equidistant  from  the  transmitter  it 
follows  that: 

(1)  If  the  attenuation  coefficient  is  not  anisotropic  (direction 

dependent)  then  the  total  signal  (H  )rp  at  each  site  due  to  the  sum  of  the 

direct  and  round-the-world  wave  will  be  identical;  The  exact  position  on  the 

interference  pattern  is  still  to  be  determined  but  the  pattern  segment  between 

the  two  receivers  will  be  symmetric  about  the  transmitter.  If  the  measured 

signal  is  the  same  we  then  can  say  that  there  is  no  directional  dependence 

for  a  even  though  we  still  have  as  yet  not  specified  what  a  really  is.  (There 

still  maybe  a  large  daytime  versus  nighttime  attenuation  rate  difference.  ) 

One  can  make  the  case  as  an  example  that  even  if  cv^,^  4  ^  may 

that  (c/v  ,  and  (c/v  ,  are  the  exact  values  which  would  make  the 

1  ph  EW  x  '  ph  WE 

signal  amplitudes  at  both  sites  equal  even  though  the  standing  wave  pattern 
is  otherwise  unsymmetrical  about  the  transmitter.  This  objection  can  be 
countered  by  changing  the  transmitter  frequency  by  a  few  Hz.  Presumably 
the  attenuation  rate  will  not  change  significantly  but  it  can  be  shown  that  the 
minima  and  maxima  positions  will  change  significantly.  If  the  signal 
amplitudes  at  this  new  frequency  are  also  the  same  then  one  must  admit 
that  there  is  no  discernible  difference  in  the  attenuation  coefficient  for  EW 
versus  WE  propagation. 
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(2)  If  the  signals  measured  are  truly  different,  then  one  must  admit 
a  significant  directional  dependence.  The  possibility  still  exists  that  one  or 
both  of  the  two  receivers  may  have  a  consistent  calibration  bias.  One  can 
circumvent  this  by  forming  the  ratios  (  |H  „/ 1 H  I  XT)  and 

(p  1,U,X,U  (P  1,U,X,1N 

(  |H  I  _ )/(  |H  I _  XT)  each  of  which  should  be  independent  of  a 

station's  calibration  bias .  As  before  if  these  ratios  are  significantly 
different  (after  the  noise  contributions  are  properly  assessed)  then  the  con¬ 
clusion  is  that  there  is  some  directional  anisotrophy.  (However  from  these 
ratios  alone  one  cannot  determine  whether  the  anisotrophy  is  associated 
with  daytime  conditions  or  with  nighttime  conditions.  )  If  on  the  other  hand 
the  ratios  are  essentially  the  same  then  the  conclusion  is  that  either  there  is 
no  anisotrophy  or  if  an  anisotrophy  exists  it  is  the  same  for  both  the  day  and 
night  conditions.  The  latter  conclusion  seems  unlikely  in  as  much  as  one 
expects  the  anisotrophy,  if  it  exists,  to  be  more  severe  at  night. 
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V. 


DATA  PRESENTATION 


The  first  point  to  be  addressed  is  the  simultaneous  calibration  of 
the  three  receiver  systems.  After  the  three  systems  were  built  and 
individually  appraised  as  being  operational,  it  was  decided  to  operate 
(at  our  Plum  Island,  Massachusetts  test  site)  all  three  systems  simultaneously 
and  compare  the  results  before  deploying  these  receivers  to  the  selected 
sites.  Unfortunately,  the  time  allocated  to  this  task  proved  insufficient. 

It  was  found  that  operating  the  three  systems  in  the  same  enclosed  area 
(a  Clark  Cortez  van)  gave  rise  to  major  interference  between  the  three 
systems  invalidating  the  results.  In  addition  several  breakdowns  on  the 
tape  recorders  occurred  and  had  to  be  repaired.  By  the  time  most  of  the 
mutual  interference  was  mitigated  and  the  recorders  repaired,  we  were  forced 
to  ship  our  equipment.  (It  should  be  noted  that  the  last  simultaneous  measure¬ 
ment  indicated  that  the  Utah  system  gave  a  1.6  dB  lower  reading  than  the  other 
remaining  system.  )  The  intent  was  to  perform  a  series  of  simultaneous 
measurements  at  Plum  Island  after  the  tests  were  finished.  At  the  beginning 
of  the  measurement  period,  it  was  realized  that  the  Utah  signal  estimates 
were  consistently  about  2  dB  too  low  in  the  40  Hz  range.  (It  should  be  pointed 
out  that  the  Utah  antenna  was  fully  buried,  the  Hawaii  antenna  half  buried 
with  volcanic  rubble  completing  the  cover  and  the  Nova  Scotia  antenna  located 
above  the  frozen  ground  in  a  specially  constructed  hut.  )  At  the  termination 
of  the  experiment  all  three  systems  were  brought  to  Plum  Island  for  a 
simultaneous  calibration.  Accordingly,  we  buried  the  Utah  antenna  to  duplicate 
conditions.  Unfortunately,  it  happened  that  the  ground  (salt)  water  level  at 
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the  first  high  tide  penetrated  the  antenna  covering  and  destroyed  the  antenna. 

Thus  we  were  left  without  a  clear  calibration  of  the  Utah  system. 

Our  response  was  to  deduce  correction  factors  by  comparing  the  Utah 

data  with  the  Nova  Scotia  data  and  with  our  previous  California  results 

(after  compensating  for  the  difference  in  distance).  Without  going  into  the 

details  of  this  comparison  we  arrived  at  correction  terms  to  be  used  in  the 

a  computation.  When  the  raw  Utah  and  Hawaii  data  are  used  in  equation  24 

|H  |Ufcah 

the  correction  terms  that  were  added  to  20  log— ^ — -  were  2.  4  dB 

^Haw 

(40  Hz  band)  and  1.0  dB  (70  Hz  band).  (It  must  be  noted  from  Fig.  6  that 
0.  5  dB  of  this  overall  correction  term  is  due  to  the  difference  in  off-axis 
corrections  for  the  Utah- Wise,  and  Hawaii- Wise,  propagation  paths.  Hence 
the  calibration  corrections  that  were  added  to  the  raw  Utah  data  are  1.9  dB 
(40  Hz  band)  and  0.  5  dB  (7  0  Hz  band).  )  It  should  be  recalled  that  the  only 
successful  simultaneous  calibration  at  Plum  Island  gave  the  Utah  system  a 
1.6  dB  lower  reading  at  45  Hz.  This  is  close  to  the  deduced  correction  value 
of  1.9  dB.  In  any  case,  the  2.  4  and  1.0  dB  correction  terms  were  the  values 
used  in  the  computation  of  the  a  estimates  which  we  have  orally  communicated 
to  the  ELF  community. 

In  October,  197  1,  during  the  following  phase  of  the  measurement 
program  we  decided  to  resolve  any  questions  concerning  our  adjustment  of 
the  Utah  data.  Accordingly,  we  simultaneously  calibrated  all  three  systems 
at  Plumlsland  (over  a  period  of  weeks).  One  of  these  systems  was  sent  to 
the  same  Utah  site  as  before  and  test  transmissions  with  the  EW  antenna 
were  scheduled.  The  results  of  these  tests  showed  with  high  confidence  the 
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TRUE  NORTH 


20  log  (cos  <!> )  =  -  1.  2  dB 


|  18  -6  -I4285~ 


Fig.  6.  'Off  axis'  pointing  angle  corrections. 
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following  to  be  true  (see  Fig.  7  for  comparison): 

(1)  The  results  with  antennas  above  or  below  ground  are  identical. 

(2)  Different  Utah  site  selections  over  a  20  mile  area  gave  the  same 
results . 

(3)  The  calibration  correction  factors  to  be  applied  to  the  earlier 
set  of  Utah  data  are  1.  9  dB  (40  Hz  range)  and  0.  8  dB  (70  Hz  range).  (This 
means  correction  terms  in  the  a  computation  of  2.4  dB  and  1.3  dB 
respectively,  when  the  off-axis  corrections  are  included.  These  numbers 
compare  extremely  well  with  the  previously  deduced  values  of  2. 4  dB  and 
1.  0  dB.  ) 

It  is  strongly  emphasized  that  the  Utah  data  in  this  report  have  been  corrected 
by  1.9  dB  and  0.  8  dB.  Off-axis  corrections  will  be  added  in  only  when  a 
computations  are  made. 

Sampling  times  and  intervals  are  chosen  to  be  the  same  for  the  Hawaii 
and  Utah  data  so  that  one  can  compare  the  signal  levels  pairwise.  The 
selection  of  a  sampling  interval  depends  on  the  signal  level  and  noise  intensity 
at  the  furthest  site  (Hawaii).  The  choice  of  the  receiver’s  integration  time  is 
somewhat  arbitrary.  It  must  be  short  enough  so  that  neither  the  amplitude 
nor  phase  of  the  received  signal  experiences  significant  change  over  the 
duration  and  long  enough  so  that  the  integrated  signal  to  noise  ratio  is 
sufficiently  large  so  that  there  is  a  reasonable  likelihood  of  detecting  true 
signal  amplitude  variations  (as  opposed  to  those  variations  caused  by  noise 
alone).  For  our  purposes,  an  integrated  signal  to  noise  ratio  of  14  dB  is  a 
desirable  lower  limit  (see  Fig.  C-2).  Dower  values  of  S/Nj-  reveal  a  rapidly 
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Fig.  7.  Comparisons  of  Utah  results  (system  calibration). 
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growing  variation  due  to  noise  alone  which  would  mask  changes  in 
signal  level.  After  experimentation  with  several  integration  times,  we  have 
selected  an  equivalent  integration  time  of  1.  42  hours,  which  we  have  used  in 
determining  all  the  estimates  in  the  following  tables  and  figures. 

The  transmission  schedules  were  selected  to  maximize  either  the 
number  of  daylight  or  nighttime  hours  over  the  entire  Nova  Scotia  to  Hawaii 
path.  No  attempt  was  made  to  determine  the  effects  introduced  by  the  day- 
night  terminator  intersecting  the  propagation  path.  One  reason  for  this 
omission  is  that  the  changes  associated  with  the  transition  from  daytime  to 
nighttime  conditions  have  a  duration  of  the  same  order  as  the  integration 
time.  Hence  accurate  signal  estimates  during  this  period  would  require  a 
significantly  shorter  integration  period,  which,  unfortunately  would  be  offset 
by  a  low  S/Nj  (with  large  variations  being  introduced  by  the  noise).  A 
generalization  that  can  be  made,  based  on  many  recording  periods  taken  all 
over  New  England,  is  that  under  either  pure  daytime  or  nighttime  path 
conditions  the  phase  change  over  a  four  to  six  hour  period  is  at  most  several 
degrees.  In  Hawaii,  the  phase  change  associated  with  day  to  night  transitions 
conditions  is  of  the  order  of  20?  Lastly,  it  appears  that  the  received  signal 
amplitude  also  is  nearly  constant  over  either  the  daytime  or  nighttime 
conditions.  There  also  exists  a  clearly  recognizable  amplitude  change  over 
the  transition  period  of  up  to  about  3  dB.  All  the  ensuing  data  should  be 
viewed  for  consistency  with  the  preceeding  generalization.  The  constancy 
of  signal  over  an  entirely  daytime  path  (or  nighttime  path)  is  the  foundation 
for  all  claims  that  we  make  that  the  average  a  and  excitation  factors  derived 
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have  a  real  physical  significance. 


Figures  8  through  15  are  typical  representatives  of  the  data  at  45  Hz  and 

7  5  Hz  at  Hawaii  and  Utah.  (The  data  at  Nova  Scotia  is  much  the  same  as  that 

2  2 

in  Utah.  )  The  vertical  scale  for  A  (or  B  )  is  proportional  to  voltage 

squared.  We  have  omitted  the  actual  voltage  scale  used  in  the  computations  and 

replaced  it  by  a  set  of  equivalent  input  field  strength  (H  )  values  in  dBH.  Note 

that  the  resulting  scales  are  very  nonlinear  in  dBH.  The  benefit  of  this 

scale  is  that  one  can  easily  see  for  example  a  *1  dB  1  variation.  The  sample 

2  2. 

values  of  and  B  are  picked  off  such  traces  and  the  biased  estimates  of 
2 

H  computed  with  the  aid  of  equation  e-8  (Appendix  E).  The  traces  for  Utah 
(D 

are  seen  to  be  relatively  constant  over  the  measurement  period.  (This  is 
consistent  with  the  S/Nj  being  generally  higher  than  24  dB.  )  The  Hawaii 
traces  exhibit  more  variability.  Much  of  this  variation  will  be  shown  to  be  due 
to  the  noise  alone;  the  S/N^  being  typically  around  12  dB  for  the  40  Hz  range 
and  16  dB  for  the  7  0  Hz  range.  It  must  be  noted  that  the  noise  intensity  often 
changes  by  as  much  as  6  dB  over  a  few  hour  period.  This  also  tends  to 
produce  traces  with  unusually  pronounced  variations  over  a  short  interval. 


The  following  table  lists  the  sample  times  chosen  for  the  two 
experiments. 
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Fig.  8.  Biased  signal  estimates 
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Biased  signal  estimates 
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Fig.  12.  Biased  signal  estimates 
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Biased  signal  estimates 
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Fig.  14.  Biased  signal  estimates 
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Fig.  15.  Biased  signal  estimates  (|H^| 


TABLE  I 


Sample  Times  (T  =  1.4 Z  hrs) 


Wisc-Utah-Hawaii  path 

For  computation 

Utah- Wisc-Nova  Scotia 

For  versus  computation 

1805  GMT 

1805 

Daytime 

1930 

1930 

2115 

2115 

2330 

2300 

0615 

0450 

Nighttime 

0748 

0615 

09  18 

0748 

1100 

0918 

These  sample  times  were  selected  to  avoid  the  effects  associated  with  the 

transition  period.  Noise  estimates  (see  Appendix  B  for  details)  are 

determined  for  the  same  sample  times  as  the  biased  signal  estimates, 

/N 

|H  L  rp.  The  unbiased  signal  estimates,  |H2|  ^  are  then  computed  as 

1  <p  !b,  T  -  1  0>u,T 

shown  in  equation  18.  Typically  the  bias  term  for  either  the  Utah  or  Nova 
Scotia  data  is  at  most  0.  1  dB.  The  noisier  Hawaii  data  exhibits  typical  bias 
terms  between  0.3  and  1.0  dB.  Figures  1 6 - Z  1  display  the  results  of  all 
the  unbiased  signal  estimates  at  all  the  sites  under  either  daytime  or  night¬ 
time  conditions.  Tables  II  through  VII  tabulate  the  same  data  along  with  the 
dates  of  transmissions. 

The  major  hypothesis  we  have  made  is  that  the  received  signal  levels 
during  a  selected  time  period  each  day  are  essentially  constant  not  only 
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Table  II.  Daytime  Unbiased  Signal  Estimates 
UTAH -NEVADA 


H  (dbH) 
<? 


T  =  1 . 42  hrs  (Equiv.  Integ.  Time) 


Tape 

Na. 

Date 

(1971) 
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(Hz) 

Sample  Times  (GMT) 

1  HJ  T 

<p  u,  T 

1805 
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2115 
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8 

Mar  8 

41 

-150.3 
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-149.8 

37 

Apr  5 

41 

-149.7 
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-150.2 
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Fig.  17.  Unbiased  signal  estimates  |  |  ^  ^ ,  Utah-Nev,  nighttime. 


50 


Table  III.  Nighttime  Unbiased  Signal  Estimates 
UTAH  -  NEVADA 


H  (db  H) 


T  =  1.42  hrs  (Equiv.  Integ.  Time) 


Tape 

Na. 

Date 

(1971) 

Freq . 
(Hz) 

Sample  Times  (GMT) 

IHJ  T 

9  u  /  T 

0450 

0615 

0748 

0918 

1100 

7 

Mar  7-8 

41 

-151.8 

-151.4 

-152.5 

-153.  1 

-152.5 

36 

Apr  4-5 

41 

-151.4 

-151.3 

-152.0 

-152.6 

-153.0 

-152.3  db 

H 

3 

Mar  15-16 

43 

-152.0 

-152.3 

-153.0 

-153.0 

-152.3 

11 

19-20 

43 

-152.9 

-152.8 

-153. 1 

-152.7 

-152. 1 

38 

Apr  5-6 

43 

-152.8 

-154.7 

-153.3 

-151.9 

-151.6 

-152.7 

1 

Mar  14-15 

45 

-153.7 

-153.4 

-152.4 

-152.0 

-152.  1 

12 

21-22 

45 

-152.5 

-153.4 

-154.0 

-153.5 

-152.0 

14 

22-23 

45 

-152.3 

-152.  1 

-155. 1 

-154.7 

-153.6 

16 

23-24 

45 

-153.7 

-153.7 

-153.4 

-152.9 

-153.  1 

34 

Apr  3-4 

45 

-152.7 

-151.7 

-152.4 

-151.8 

-150.4 

-152.8 

5 

Mar  16-17 

47 

-151.7 

-151.4 

-152.7 

-152.0 

-152.3 

18 

24-25 

47 

-151.8 

-152.0 

-152.6 

-153.1 

-152.5 

20 

25-26 

47 

-152.3 

-152.3 

-152.6 

-152.6 

-152.3 

-152.4 

9 

Mar  18-19 

49 

-151.3 

-151.3 

-151.7 

-151.6 

-151.3 

40 

Apr  6-7 

49 

-151.2 

-151.6 

-151.5 

-151.7 

-151.4 

-151.5 

26 

Mar  30-31 

73 

-149.6 

-149.7 

-149.9 

-148.9 

-148.4 

28 

Mar  31  - 

Apr  1 

73 

-148.3 

-148.5 

-149.9 

-148.7 

-148.  1 

-149.0 

22 

Mar  28-29 

75 

-148.2 

-147.9 

-148.4 

-149. 1 

-149.2 

24 

Mar  29-30 

75 

-148.5 

-148.2 

-149.2 

-150.0 

-148.4 

-148.9 

30 

Apr  1-2 

77 

-147.8 

-147.5 

-148.0 

-148.6 

-148.4 

32 

-  1 

2-3 

77 

-149.  1 

-149.0 

-149.7 

-149.2 

-148.6 

-148.6 
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SAMPLE  TIMES  (GMT) 
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Table  IV.  Daytime  Unbiased  Signal  Estimates 
HAWAII 


H  (dbH) 


T  =  1 . 42  hrs  (Equiv.  Integ,  Time) 


Tape 

Na. 

Date 

(1971) 

Freq. 

(Hz) 

Sample  Times  (GMT) 

ITTl  _ 

t  u,  T 

1805 

1930 

2115 

2300 

2330 

8 

Mar  8 

41 

-158.4 

-161.0 

-160.6 

-159.5 

37 

Apr  5 

41 

-157.2 

-158.6 

-158.6 

-158.8 

-159.2 

4 

Mar  16 

43 

-162.3 

-162.7 

-161.2 

-158.6 

39 

Apr  6 

43 

-161.4 

-158.1 

-160.4 

-162.6 

-160.5 

2 

Mar  15 

45 

-158.9 

-157.3 

-159.9 

-158.  1 

13 

22 

45 

-162.2 

-158.3 

-157.9 

-159.0 

15 

23 

45 

-159.  1 

-158.6 

-159.2 

-158.4 

17 

24 

45 

-160.  1 

-159.6 

-159.6 

-158.8 

35 

Apr  4 

45 

-159.2 

-157.6 

-157.7 

-159.4 

-158.8 

6 

Mar  1 7 

47 

-158.4 

-158.0 

-158.8 

-157.8 

19 

25 

47 

-158.6 

-157.7 

-158.8 

21 

26 

47 

-156.2 

-155.8 

-156.7 

-157.5 

-157.8 

10 

Mar  19 

49 

-158.8 

-157.5 

-158.2 

-156.7 

41 

Apr  7 

49 

-157.6 

-155.0 

-156.8 

-155.5 

-157.2 

27 

Mar  31 

73 

-158.3 

-158.2 

-156.6 

-157.4 

29 

Apr  1 

73 

-157.9 

-159.  1 

-159.0 

-157.9 

-157.8 

23 

Mar  29 

75 

-157.5 

-156.4 

-157.3 

-157.7 

25 

30 

75 

-157.3 

-157.5 

-158.8 

-158.9 

-157.6 

31 

Apr  2 

77 

-157.7 

-156.3 

-155.8 

-156.9 

33 

3 

77 

-156.4 

-156.4 

-156.5 

-156.5 

-156.6 

53 
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Fig.  19.  Unbiased  signal  estimates  |H  ^  ,  Hawaii,  nighttime. 
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Table  V.  Nighttime  Unbiased  Signal  Estimates 
HAWAII 


H  (dbH) 
9 


T  =  1.42  hrs  (Equiv.  Integ.  Time) 


Tape 

Na. 

Date 

(1971) 

Freq. 

(Hz) 

Sample  Times  (GMT) 

lH~l  T 

t  u,T 

0450 

0615 

0748 

0918 

1100 

7 

Mar  7-8 

41 

-162.0 

-163.5 

-161.5 

-161.6 

36 

Apr  4-5 

41 

-159.0 

-159.0 

-163. 1 

-161.7 

-161.2 

3 

Mar  15-16 

43 

11 

19-20 

43 

-158.4 

-158.0 

-157.7 

-160.9 

38 

Apr  5-6 

43 

-160.7 

-161.4 

-159.0 

-159.7 

-159.3 

1 

Mar  14-15 

45 

-155.8 

-159. 1 

-159.5 

-158.6 

12 

21-22 

45 

-157.2 

-158.9 

-161.1 

-160.3 

14 

22-23 

45 

-159.6 

-153. 1 

-162.6 

-161.4 

(-159.4)  - 

-  16  pt  ave 

16 

23-24 

45 

-159.5 

-161.9 

-160.6 

-161.4 

34 

Apr  3-4 

45 

-158.0 

-157.  1 

-156.  1 

-155.3 

-158.7  - 

-  20  pt  ave 

5 

Mar  16-17 

47 

-159.5 

-160.7 

-164.2 

-163.4 

18 

24-25 

47 

-158.8 

-161.6 

-161.0 

-161.0 

20 

25-26 

47 

-159.5 

-160.2 

-161.1 

-162.5 

-160.7 

9 

Mar  18-19 

49 

-161.0 

-160.2 

-163.8 

-164.2 

40 

Apr  6-7 

49 

-158.9 

-159.7 

-159.4 

-161.0 

-160.8 

26 

Mar  30-31 

73 

-157.2 

-158.4 

-158.3 

-158.3 

28 

Mar  31- 

Apr  1 

73 

-157.7 

-158.  1 

-158.2 

-157.6 

-157.8 

22 

Mar  28-29 

75 

-158.9 

-158.  1 

-157.2 

-158. 1 

24 

Mar  29-30 

75 

-154.9 

-155.2 

-158.4 

-157.3 

-157.0 

30 

Apr  1  -2 

77 

-157.0 

-155.9 

-156.4 

-155.9 

32 

2-3 

77 

-156.0 

-158.2 

-156.6 

-155.7 

-156.4 
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Table  VI.  Daytime  Unbiased  Signal  Estimates 
NOVA  SCOTIA 


H  (dbH) 
9 


T-  1.42  hrs  (Equiv.  Integ.  Time) 


Tape 

Na. 

Date 

(1971) 

Freq. 

(Hz) 

Sample  Times  (GMT) 

1805 

1930 

21 15 

2300 

2330 

lH  (fj  T 
<P  1  “,T 

8 

41 

27 

Apr  5 

41 

-149.4 

-149.8 

-149.8 

-150.  1 

-149.8 

4 

Mar  16 

43 

-150.0 

-149.2 

-149.0 

-149.9 

39 

Apr  6 

43 

-148.6 

-148.3 

-148.6 

-148.5 

-149.0 

2 

Mar  15 

45 

-149.  1 

-149.0 

-149.2 

-149.2 

13 

22 

45 

-148.2* 

-147.8* 

-148.8 

-149.  1 

15 

23 

45 

-147.7* 

-147.6* 

-147.7 

-149.  1 

17 

24 

45 

35 

Apr  4 

45 

-147.6* 

-148.4 

-148.5 

-148.5 

-148.4 

6 

Mar  17 

47 

19 

25 

47 

-148. 1 

-148.  1 

21 

26 

47 

-148.4 

-148.5 

-149.  1 

-149.0 

-148.5 

10 

Mar  19 

49 

-149.0 

-148.7 

-149.0 

-149.0 

41 

Apr  7 

49 

-147.  1 

-147.1 

-147.5 

-148.4 

-148.2 

27 

Mar  31 

73 

t 

29 

Apr  1 

73 

_ t 

23 

Mar  29 

75 

-145.9 

-145.6 

-146.3 

-146.9 

25 

30 

75 

-146.2 

-146.7 

t 

-146.3 

31 

Apr  2 

77 

t 

33 

3 

77 

t 

-146.2 

-145.9 

-146.0 

*  Naisy  period, 
t  Tape  recorder  problems. 
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Fig.  21.  Unbiased  signal  estimates  |H  Nova  Scotia,  nighttime. 
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Table  VII.  Nighttime  Unbiased  Signal  Estimates 
NOVA  SCOTIA 


H  (db  H) 


T=  1.42  hrs  (Equiv.  Integ.  Time) 


Tape 

Na. 

Date 

(1971) 

Freq. 

(Hz) 

Sample  Times  (GMT) 

0450 

0615 

0748 

0918 

1100 

|H  (f  )| 

<p  \  u,\ 

7 

Mar  7-8 

41 

36 

Apr  4-5 

41 

-152.8 

-152.9 

-152.2 

-151.0 

-152.2 

3 

Mar  15-16 

43 

-151.9 

-152.  1 

-152.4 

-151.7 

11 

19-20 

43 

-154.7 

-154.7 

-154.9 

-153.3 

38 

Apr  5-6 

43 

-152.3 

-152.8 

-151.5 

-150.4 

-152.6 

1 

Mar  14-15 

45 

-151.0 

-151.0 

12 

21-22 

45 

-152.9 

-152.2 

-151.4 

-151. 1 

14 

22-23 

45 

-152.2 

-153.2 

-150.8 

-151.0 

16 

23-24 

45 

-153.5 

-153.5 

34 

Apr  3-4 

45 

-152.  1 

-150.  1 

-151.0 

-149.8 

-151.6 

5 

Mar  16-17 

47 

-151.4 

-151.1 

-150.8 

-151.5 

18 

24-25 

47 

-153.0 

-152.8 

-153.9 

-153.7 

20 

25-26 

47 

-151.8 

-153.7 

-153.6 

-153.1 

-152.5 

9 

Mar  18-19 

49 

-151.5 

-151.5 

-151.0 

-151.3 

40 

Apr  6-7 

49 

-151.8 

-151.1 

-150.5 

-150.6 

-151.2 

26 

Mar  30-31 

73 

-149.7 

-148.8 

-149.2 

-148.9 

-149.  1 

28 

Mar  31  - 

Apr  1 

73 

22 

Mar  28-29 

75 

-147.9 

-148,0 

-148.9 

-149.8 

24 

Mar  29-30 

75 

-147.9 

-147.6 

-148.6 

-148.3 

-148.3 

30 

Apr  1-2 

77 

32 

2-3 

77 

-147.4 

-148.2 

-148.2 

-147.2 

-147.7 

59 


over  the  period  but  also  over  similar  time  intervals  on  succeeding  days. 

If  the  premise  is  correct  then  the  &  one  obtains  by  suitably  averaging  all 
the  data  has  a  real  physical  meaning.  If  the  premise  is  not  quite  correct 
then  these  computed  values  of  q,  are  Average  1  values  for  the  attenuation 
rates.  The  true  attenuation  rate  for  a  particular  time  interval  could  then 
be  a  few  tenths  of  a  db/Mm  different.  With  the  transmitted  power  levels 
currently  available  with  WTF,  it  is  clear  that  we  will  always  require  a 
large  number  of  samples  (with  long  integration  times  per  sample)  in  order 
to  calculate  an  average  value  for  #  (with  reasonable  error  bounds).  We 
can,  however,  look  at  the  variability  of  received  signal  estimates  at  the  more 
remote  receiver  sites.  Examination  of  these  estimates  along  with  their 
individual  confidence  limits  may  provide  assurance  that  the  hypothesis  made 
is  reasonable.  In  order  to  compute  the  confidence  limits  on  each  signal 
estimate,  the  corresponding  noise  estimates  are  also  computed  at  the  same 
sample  times.  As  explained  in  Appendix  B,  about  2400  independent  samples 
(over  a  one  hour  period)  are  taken  in  order  to  compute  each  noise  intensity 
value.  The  noise  estimates  are  expected  to  be  accurate  within  ±  .15  db 
provided  the  noise  is  reasonably  stationary  over  the  one  hour  period.  As 
typical  examples,  tables  VIII  and  IX  list  the  values  of  noise  intensity  and  rms 
signal  to  integrated  noise  ratios  for  Hawaii.  For  each  signal  to  integrated 


)  confidence  intervals  are 


constructed  (see  appendix  C  for  details).  Figures  22  through  25  show  the 


unbiased  signal  estimates  and  the  80<$  confidence  intervals  for  the  Utah 


and  Hawaii  data  at  45,  73,  7  5  and  77  Hz.  Some  important  conclusions  can 
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Table  VIII.  Doytime  Noise  Intensities  ond  Integroted  Signal  to  Noise  Ratios 


HAWAII 


dbH/VHl  Sb/N(  N|=Hn/t  T  =  1 . 42  hrs 


Tape 

Na. 

Dote 

(1971) 

Freq. 

(Hz) 

Sample  Times  (GMT) 

Sample  Times  (GMT) 

1805 

1930 

2115 

2330 

1805 

1930 

2115 

2330 

8 

Mar  8 

41 

-133.3 

-133.4 

-137.0 

-137.3 

12.3 

10.0 

13.7 

15.  1 

37 

Apr  5 

41 

-133.0 

-133.0 

-134.3 

-134.3 

13.1 

11.8 

13.0 

12.8 

4 

Mar  16 

43 

-134.4 

-134.4 

-134.8 

-135.2 

9.7 

9.3 

11.1 

13.9 

39 

Apr  6 

43 

-134.0 

-134.0 

-133.9 

-133.9 

10.  1 

13.2 

11.0 

9.0 

2 

Mar  15 

45 

-133.6 

-133.6 

-135.3 

-135.3 

12.  1 

13.6 

12.7 

14.5 

13 

22 

45 

-132.5 

-132.5 

-134.2 

-135.9 

8.  1 

11.6 

13.6 

14.2 

15 

23 

45 

-133.9 

-133.9 

-134.3 

-136.5 

12.2 

12.6 

12.5 

15.3 

17 

24 

45 

-134.4 

-134.4 

-136.  1 

-138.6 

11.7 

12.3 

13.8 

17.0 

35 

Apr  4 

45 

-133.  1 

-133.6 

-134.3 

-134.3 

11.3 

13.3 

13.9 

12.3 

6 

Mar  17 

47 

-135.5 

-135.5 

-136.7 

-136.7 

14.4 

14.7 

15.1 

16.  1 

19 

25 

47 

-132.7 

-132.7 

-132.7 

-133.2 

11.5 

12. 1 

11.8 

21 

26 

47 

-132.0 

-131.3 

-132.1 

-132.8 

13.  1 

12.8 

12.7 

12.6 

10 

19 

49 

-136.0 

-136.0 

-136.5 

-137.0 

14.5 

15.7 

15.5 

17.5 

41 

Apr  7 

49 

-132.0 

-132.0 

-131.5 

-131.5 

11.8 

14.3 

12.  1 

13.3 

27 

Mar  31 

73 

-138.7 

-138.7 

-141.5 

-142.4 

17.6 

17.7 

22.0 

22.  1 

29 

Apr  1 

73 

-139.  1 

-139.  1 

-140.8 

-140.8 

18.4 

17.2 

19.0 

20.0 

23 

Mar  29 

75 

-139.  1 

-139.  1 

-141.8 

-141.8 

18.8 

19.8 

21.6 

21.2 

25 

30 

75 

-135.5 

-135.5 

-137.9 

-137.9 

15.4 

15.  1 

16.3 

16.2 

31 

Apr  2 

77 

-136.4 

-136.4 

-136.2 

-136.2 

15.9 

17.2 

17.6 

16.5 

33 

3 

77 

-136.4 

_ 

-136.3 

-139.  1 

-139.1 

17.2 

17.0 

19.7 

19.7 
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Table  IX.  Nighttime  Naise  Intensities  and  Integrated  Signal  ta  Naise  Ratios 


HAWAII 


H^|  dbH A/Hi  S/Nj  Nj  =  H^/T  T=1.42hrs 


Tape 

Na. 

Date 

(1971) 

Freq. 

(Hz) 

Sample  Times  (GMT) 

Sample  Times  (GMT) 

0615 

0748 

0918 

1100 

0615 

0748 

0918 

1100 

7 

Mar  7-8 

41 

-135.6 

-135.6 

-135.6 

-134.3 

11.1 

9.7 

11.5 

10.2 

36 

Apr  4-5 

41 

-134.7 

-134.7 

-134.7 

-134.3 

13.0 

13.0 

9.3 

10.  1 

3 

Mar  15-16 

43 

11 

19-20 

43 

-134.4 

-133.4 

-134.4 

-135.5 

13.3 

12.7 

14.0 

12.0 

38 

Apr  5-6 

43 

-133.6 

-133.6 

-133.0 

-132.9 

10.4 

9.8 

11.4 

10.7 

1 

Mar  14-15 

45 

-133.1 

-133.1 

-134.5 

-134.5 

14.5 

11.4 

12.4 

13.2 

12 

21-22 

45 

-133.6 

-134.8 

-136.4 

-136.4 

13.7 

13.2 

12.6 

13.4 

14 

22-23 

45 

-132.9 

-132.9 

-133.3 

-133.3 

10.8 

12.2 

9.0 

9.5 

16 

23-24 

45 

-135.6 

-135.8 

-134.5 

-134.5 

13.4 

11.3 

.  11,3 

10.6 

34 

Apr  3-4 

45 

-133.4 

-133.4 

-134.  1 

-134.  1 

12.7 

13.6 

15.2 

16.0 

5 

Mar  16-17 

47 

-136.6 

-136.6 

-135.2 

-135.2 

14.4 

13.2 

8.7 

9.4 

18 

24-25 

47 

-136.7 

-136.7 

-136.7 

-135. 1 

15.  1 

12.5 

13.0 

11,5 

20 

25-26 

47 

-134.0 

-134.0 

-133.0 

-131.0 

11.9 

11.2 

9.5 

7.4 

9 

Mar  18-19 

49 

-134.8 

-134.8 

-135.2 

-135.4 

11.2 

12.0 

9.  1 

8.9 

40 

Apr  6-7 

49 

-134.6 

-134.6 

-134.8 

-135.0 

13.0 

12.3 

12.7 

11.4 

26 

Mar  30-31 

73 

-137.0 

-138.3 

-138.4 

-138.5 

17.0 

17.  1 

17.3 

17.4 

28 

Mar  31  - 

Apr  1 

73 

-139.0 

-138.1 

-138.  1 

-138.9 

18.5 

17.2 

18.  1 

18.5 

22 

Mar  28-29 

75 

-136.4 

-136.8 

-136.8 

-137.4 

14.7 

15.9 

16.8 

16.5 

24 

29-30 

75 

-137.3 

-136.5 

-136.9 

-137.3 

19.5 

18.5 

15.7 

17.2 

30 

Apr  1  -2 

77 

-134.0 

-135.4 

-136.0 

-136.4 

14.3 

16.7 

16.8 

17.7 

32 

i 

2-3 

77 

-133.9 

-133.9 

-134.5 

-135.5 

15.  1 

13.0 

15.  1 

17.0 

62 


{ d  B  H  ) 


Fig.  22.  Unbiased  signal  estimates 
confidence  intervals,  daytime. 
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Fig.  24.  Unbiased  signal  estimates 
confidence  intervals,  nighttime. 
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be  drawn  from  these  results.  For  the  daytime  data  (figures  22  and  23) 
the  confidence  intervals  for  the  Utah  and  Hawaii  data  have  a  significantly 
large  common  overlap  interval.  This  fact  very  strongly  suggests  that  the 
daytime  signal  strengths  from  sample  to  sample  (and  from  day  to  day) 
are  nearly  constant.  Implicit  in  this  statement  is  that  both  the  attenuation 
rate  and  excitation  factor  (in  both  the  40  Hz  and  70  Hz  ranges)  are  to  a 
high  degree  constant  quantities.  It  is  then  felt  that  the  averaging  process 
to  be  undertaken  will  give  us  average  attenuation  rates  and  excitation 
factors  close  to  the  true  values.  The  nighttime  data  (figures  24  and  2  5) 
exhibit  significantly  more  variability  even  though  the  integrated  signal  to 
noise  ratios  are  not  that  different.  It  would  be  difficult  here  to  argue 
that  the  signal  level  is  in  fact  a  constant.  We  maintain  however  that  the 
signal  level  changes  are  due  primarily  to  changes  in  the  excitation  factor. 
This  is  quite  clearly  seen  in  examination  of  figure  2  5.  It  is  noted  that  the 
signal  variations  between  Utah  and  Hawaii  are  highly  correlated.  This 
is  consistent  with  an  excitation  factor  change  since  such  a  change  effects 
both  sites  equally  and  this  appears  to  be  the  case.  Extremely  large  changes 
in  &  would  be  required  to  fit  the  data  at  Utah  if  one  attributed  signal  level 
variations  as  being  due  primarily  to  #  changes.  Moreover,  this  would  be 
reflected  in  even  larger  signal  variations  in  Hawaii.  This  appears  not  to 
be  true.  Our  position  then  is  that  the  signal  variations  at  night  are  mostly 
attributable  to  variations  in  the  excitation  factor.  Since  the  excitation 
factor  cancels  out  when  we  combine  the  data  for  the  two  sites,  the  average 
nighttime  attenuation  rate  computed  is  a  reasonable  approximation  to  the 
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actual  nighttime  attenuation  rates. 


The  size  of  the  confidence  intervals  in  figures  22  through  25  il¬ 
lustrates  an  extremely  important  point.  That  is,  any  attempt  to  compute 
an  instantaneous  'real'  value  of  a  from  any  single  pair  of  signal  estimates 
is  going  to  be  totally  useless  because  of  the  large  uncertainty  introduced 
by  noise.  Moreover,  in  order  to  measure  a  relatively  noise  free 
instantaneous  'real [  values  of  #  (with  integration  times  of  the  order  of 
30  minutes  or  less),  one  would  need  to  increase  the  radiated  power  by 
at  least  20  db.  This  is  clearly  very  expensive  and  hence  for  the  foreseeable 
future,  one  is  going  to  be  forced  to  rely  on  average  estimates  of  #  as  well 
as  excitation  factor. 

The  variation  in  the  nighttime  excitation  factor  also  indicates  the 
absolute  necessity  of  making  multi  site  measurements  if  one  wants  to 
deduce  reasonable  estimates  of  q,.  Single  site  measurements  have  their 
value  in  the  estimation  of  signal  strengths  for  that  particular  locations. 

One  could  over  a  long  period  determine  the  statistics  of  the  signal  variation^ 
however,  relating  these  statistics  to  the  propagation  parameters  would  be 
exceedingly  difficult. 

During  October  1971,  Lincoln  Laboratory  performed  another  series 
of  ELF  measurements  with  sites  in  Utah,  Greenland  and  Norway  (the  re¬ 
sults  of  the  experiment  are  still  undergoing  processing).  The  preliminary 
conclusions  are  essentially  identical  to  those  discussed  in  this  report. 

Figures  26  and  27  exhibit  the  final  signal  averages  (and  confidence 
intervals)  as  a  function  of  site  location  and  frequency.  The  method  of 
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Fig.  26.  Average  field  strengths,  |H^(f.)|  and  |H^|  (daytime), 
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construction  of  confidence  curves  and  hence  the  determination  of  confidence 


intervals  for  the  averages  (  |H^  (L)  |T)  is  explained  fully  in  appendix  D. 
Figure  D-I  shows  one  such  set  of  curves  for  the  day  and  night  averages  of 
the  Hawaii  data  at  45  Hz.  The  averages  (  |H^  (f.)  |T)  were  computed  from 
the  unbiased  signal  estimates  and  noise  estimates  through  equation  21. 
These  averages  are  also  tabulated  in  Tables  II  through  VII. 

In  order  to  reduce  the  effects  of  the  standing  wave  on  the  attenuation 
rate  computation,  it  has  been  decided  to  perform  a  second  average  (over 


on 


frequency)  on  the  earlier  single  frequency  averages  [H^  (f.)  ,p.  Equati 

22  shows  the  appropriate  weighting.  The  result  IH^  (f  )  I  is  computed 

for  each  site  for  both  daytime  and  nighttime  conditions. 


TABLE  X 


|Hjp  (fm)  1  u>  Estimates 

(dbH) 

40  Hz  Band 

f  =  45  Hz 
m 

70  Hz  Band  f  =75  Hz 
m 

Day 

Night 

Day 

Night 

Nova  Scotia 

-148. 8 

-152. 0 

-146.1 

-148. 3 

Utah 

-149. 3 

-152.3 

-146. 6 

-148. 8 

Hawaii 

-158. 6 

-160. 0 

-157. 3 

-157.0 

The  80$  confidence  intervals  for  the  values  of  |Hm  (f  )  I  were  computed 

r  I  H-'  m  lU,  $ 

as  explained  in  Appendix  D.  The  confidence  interval  on  the  Utah  values  are 
less  than  ±  0.1  db  but  are  not  displayed  since  other  sources  of  error  (i.e.  , 
other  than  noise)  would  dominate.  It  is  noteworthy  that  the  confidence 
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interval  on  the  Hawaii  averages  is  of  the  order  of  ±  0.3  db/Mm,  This 
would  imply  that  the  values  of  a  could  be  computed  to  less  than  ±  0.1 
db/ Mm  provided  the  signal  strengths  are  in  fact  true  constant  levels  and 
provided  that  other  sources  of  calibration  error  are  not  introduced. 

The  confidence  intervals  for  the  Nova  Scotia  data  are  approximately 
the  same  as  those  on  Utah.  However,  throughout  the  experimental  period, 
the  tape  recorder  exhibited  speed  variations  which  has  affected  many  of 
the  signal  estimates.  This  effect  is  generally  small  but  probably  of  the 
same  order  as  the  confidence  interval.  We  maintain  that  this  speed  con¬ 
trol  error  on  the  Nova  Scotia  tape  recorder  accounts  for  some  of  the 
variation  in  differences  in  signal  levels  as  measured  in  Utah  and  Nova  Scotia. 
From  figure  6  one  observes  that  if  and  were  equal,  then  due  to 

the  off  axis  correction, the  Nova  Scotia  averages  should  be  0.6  db  higher 
than  those  in  Utah.  Examination  of  figures  26  and  27  show  that  this  is 
very  likely.  In  any  case  it  would  be  consistent  with  the  data  to  state  that 
the  attenuation  rate  difference  (i.e.  ,  Ic^jr  "  |  )  less  than  0.2 

db/Mm  for  both  the  40  Hz  and  the  70  Hz  range.  (The  final  averages 
(over  frequency)  in  figures  26  and  27  show  consistent  differences  of  about 
0.  5  db  at  both  the  40  Hz  and  70  Hz  bands  for  day  and  night.  We  maintain 
that  if  differences  in  -  &EW  I  they  are  so  insignificant  that 

we  are  unable  to  measure  them  with  the  transmitter  power  available. 

Equation  1  exhibits  a  direct  f  dependence  in  addition  to  the  weaker 

implicit  dependence  of  f  in  a',  c/v  ^  ,  h  and  (T  ^  .  Accordingly,  due  to  the 

49 

direct  dependence  alone,  one  would  expect  a  1.  5  db  (20  log  yy-  )  increase 
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across  the  40  Hz  band  and  a  0.  5  db  increase  across  the  70  Hz  range.  Ex¬ 
amination  of  figures  26  and  27  show  such  increases  with  a  few  exceptions. 
The  daytime  Hawaii  results  in  the  40  Hz  range  would  show  the  1.  5  db 
increase  except  for  the  lower  than  expected  average  at  43  Hz.  Figure  18 
shows  the  collection  of  unusually  low  values  of  signal  strength  which  re¬ 
sulted  in  this  average.  The  nighttime  averaged  results  (Fig.  27)  in  the 
40  Hz  range  are  an  extremely  poor  fit  to  a  1.5  dB  increase  over  this  band. 
Even  the  inclusion  of  the  effects  of  round-the-world  interference  would  not 
make  the  fit  appreciably  better.  It  can  be  shown  using  the  attenuation  rates 
listed  in  Table  XI  that  the  maximum  variation  in  signal  level  at  Hawaii 
would  only  be  ±  0.  3  dB.  This  is  not  enough  to  provide  a  reasonable  fit  to 
the  data.  It  must  be  acknowledged  simply  that  the  nighttime  signals  show 
more  variability. 


73 


VI.  CONCLUSIONS 


From  Table  X  and  equation  2  4  one  can  finally  compute  the 
attenuation  rates.  The  resulting  values  are  shown  in  Table  XI.  The 
significant  point  is  that  the  daytime  attenuation  rates  are  higher  than  earlier 
theoretical  predictions  indicated.  The  question  of  error  bars  on  the  average 
values  of  a  is  answered  in  the  following  way.  From  the  preceding  section 
it  was  shown  that  the  error  bars  (confidence  intervals)  on  a  due  to  the 
atmospheric  noise  alone  are  less  than  ±  0.  1  dB/Mm  (provided  the  signal  level 
is  really  a  constant  value).  The  error  budget  due  to  system  calibration 
errors  (see  Appendix  I)  shows  maximum  error  bounds  on  the  signal  estimates 
of  ±  0.2  dB  per  site.  Hence  the  error  bounds  on  a  due  to  system  calibra¬ 
tion  errors  are  again  less  than  ±  0.  1  dB/Mm.  The  total  confidence  interval 
for  the  a  measurements  is  then  about  ±  0.  15  dB/Mm.  We  maintain  that  the 
average  daytime  attenuation  rates  are  in  fact  extremely  close  to  the 
instantaneous  values  of  o'.  This  statement  is  based  on  the  relatively  constant 
signal  levels  during  the  day.  We  note  that  our  samples  come  from  a  one 
month  period  when  the  solar  conditions  were  reasonably  quiet.  The  nighttime 
signal  levels  are  noticeably  more  variable.  A  good  deal  of  this  variation 
is  due  to  the  changes  in  excitation  function  (that  is  variations  in  the 
equivalent  height  of  the  ionosphere).  The  variation  on  the  instantaneous 
value  of  nighttime  a  is  at  the  present  not  measurable.  Our  current 
speculation  is  that  the  variation  is  perhaps  of  the  order  of  ±  0.2  dB/Mm. 

The  ratio  of  the  daytime  to  nighttime  excitation  factors  is  computed 
from  equation  29.  The  results  are  tabulated  in  Table  XI.  If  one  uses  the 
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Table  XI.  Attenuation  Rates,  Excitation  Factors,  Ionospheric  Heights 


o  (db/Mm) 

0) 

HP 

(3) 

h  (km) 

(4) 

E  = - - - 

HN  (2) 

hA(c/vph)  (5) 

45  Hz 

Day 

1.  16 

1.53 

(3.68  db) 

1.26 

69.2 

-58.60  db 

Night 

0.81 

1.07 

113.0 

-62.28 

75  Hz 

Day 

1.47 

1.52 

(3.64  db) 

1.28 

65.7 

-59.72 

Night 

0.92 

1.20 

103. 

-63.36 

20  log  0.5 

I'Ll  H 

(1)  From  equotion  24  where  a  =  —  .99  + -  r  - 

D  4. 56 


.5  is  the  off  axis  correction 


(2)  From  equotion  29 

(3)  From  U5C/NLL 

(4)  Depends  on  volidity  of  a  and 

(5)  From  equation  23  after  o  determined 


c/v  ,  estimates 
ph 
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measured  values  of 


<j  ^  and  some  postulated  values  of  c/v  the  values  of 
the  equivalent  ionospheric  heights  can  be  determined.  For  the  daytime  in 
the  40  Hz  band,  the  equivalent  ionosphere  is  69  km.  This  is  consistent  with 
the  range  of  previous  theoretical  estimates  for  h.  The  slightly  lower  value 
of  65.7  km  for  the  70  Hz  band  is  also  reasonable  (the  lower  value  probably 
being  due  to  a  poor  choice  in  c/v  j^)*  The  nighttime  ionospheric  heights 
are  113  and  103  km  for  the  40  and  70  Hz  bands  respectively.  A  priori,  one 
would  have  expected  values  closer  to  90  km.  The  nighttime  excitation 
factor  is  then  seen  to  be  respectively  2.  0  and  1. 2  dB  worse  than  expected. 

Finally,  with  the  most  accurate  measurements  that  we  are  at  present 
capable  of  performing,  no  apparent  difference  in  -  c^p  |  seem 

Certainly  this  difference  is  less  than  ±  .2  dB/Mm. 

Table  XII  is  constructed  with  the  aid  of  the  measured  values  of  a  and 
^d/^N'  The  totals  agree  well  with  the  measured  averages  in  Table  X  (there 
is,  of  course,  some  commonality  in  that  Table  X  is  used  to  determine  a  and 
p  and  then  these  are  used  in  Table  XII).  However,  the  consistency  of  both 
tables  leads  one  to  believe  that  equation  1  is  in  fact  a  good  formula  to  des¬ 
cribe  the  propagation  at  these  frequencies  and  distances. 
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Table  XII.  Signal  Estimates 
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X 

r>. 
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APPENDIX  A 


Equipment  Design  Considerations 

The  equipment  used  in  the  propagation  experiments  can  logically  be 
separated  into  the  field  recording  equipment  and  the  laboratory  playback 
facility.  What  follows  is  a  description  of  the  field  recording  equipment  and 
a  discussion  of  the  motivation  behind  the  design. 

The  field  equipment  consists  of  three  units:  the  antenna,  two  racks  of 
electronic  equipment,  and  a  D.  C.  motor  generator  for  prime  power  genera¬ 
tion.  The  antenna,  a  5  foot  loop,  is  either  supported  vertically  in  a  hut  or 
buried  in  the  ground  and  is  separated  approximately  ZOO  feet  from  the  elec¬ 
tronics  racks.  Approximately  250  feet  separate  the  motor  generator  set 
and  the  antenna.  All  site  locations  have  been  selected  to  be  as  far  as  possi¬ 
ble  from  power  lines,  typically  this  distance  is  in  the  order  of  3  -  5  miles. 
Whenever  possible,  the  site  is  selected  only  after  preliminary  exploration 
with  a  M  sniffer"  loop  to  determine  the  extent  of  power  line  contamination. 
Interference  at  a  site  is  influenced  by  the  geometry  of  the  power  lines  and 
pointing  angle  of  the  antenna  as  well  as  separation. 

The  pre-amplifier,  filters,  tape  recorder  and  other  items  are  mounted 
in  two  19!I  racks  26  inches  high.  On  site  these  racks  are  mounted  in  a  panel 
type  van.  Under  cold  weather  conditions  the  equipment  is  placed  in  an 
environmental  chamber  equipped  with  thermostatic  control  to  maintain  the 
temperature  in  the  60  -  70°F  range.  A  26  inch  rack  height  limitation  was 
selected  for  two  reasons:  (1)  these  racks  will  fit  into  the  rear  of  a  station 
wagon,  if  necessary,  and  (2)  they  are  easily  air  transportable. 
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The  D.  C.  motor  generator  is  a  commercially  available  combination 
which  has  been  integrated  with  a  larger  gas  tank  into  an  environmentally 
protected  stand  alone  unit. 

A  discussion  of  the  field  recording  equipment  follows. 

Antenna  Design 

The  choice  of  a  particular  antenna  configuration  was  influenced  by 
several  considerations.  Perhaps  the  earliest  choice  that  had  to  be  made 
was  that  between  a  "cored”  antenna  or  an  air  core.  It  is  obvious  that  a  loop 
with  a  magnetic  core  can  be  made  smaller  for  a  given  turns  area  product. 

It  has  been  our  experience,  however,  that  the  permeability  of  a  magnetic 
core  loop  is  subject  to  change  due  to  various  effects.  Depending  on  the  type 
of  core  material  and  the  physical  configuration  sensitive  antennas  such  as 
these  may  have  their  permeability  altered  by  the  earth1  s  magnetic  field. 

This  effect  can  be  compensated  for;  however,  it  would  introduce  another 
possible  source  of  error.  Additionally,  we  have  observed  at  least  one 
core  type  antenna  whose  NA  product  has  changed  by  a  factor  of  two  over  a 
period  of  two  years.  Whether  this  was  due  to  shock  while  being  shipped, 
abnormal  current  surges,  or  mechanical  failure  is  not  known,  but  it  does 
re-inforce  the  conviction  that  a  stable  air-core  loop  is  best  suited  for  pre¬ 
cision  measurements. 

The  open  circuit  voltage  generated  in  an  air  core  loop  antenna  due  to 
a  magnetic  field  is  given  by: 

V  =  |h  U  2tt  f  N  A  a -1 

oc  1  <p 1  o 
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H  -  magnetic  field,  amp/meter 

Ijl^  -  permeability  of  free  space,  4n  x  10  henry/meter 

f  -  frequency  of  applied  field  assumed  sinusoidal  (hertz) 

N  -  number  of  turns  on  antenna 

2 

A  -  area  of  loop  (meters) 

From  a  measurement  point  of  view  the  two  antenna  parameters  which 
are  critical  to  accurate  field  estimates  are  the  number  of  turns  and  the  area 
of  the  loop.  These  two  parameters  are  often  lumped  together  and  referred 
to  as  the  turns-area  product  (NA).  The  choice  of  a  large  diameter  loop 
with  fewer  turns  over  a  smaller  diameter  with  many  turns  was  made  on  the 
basis  that  the  larger  diameter  would  allow  a  more  precise  estimate  of  the 
turns-area  product  thus  introducing  less  error  into  the  measurements. 

The  antenna  is  constructed  with  471  turns  of  no.  14  double  formvar 
insulated  copper  wire  wound  within  a  loop  of  5.  862  cm.  OD  Poly  Vinyl 
Chloride  (PVC)  tubing  formed  into  a  torus  with  inside  diameter  of  1.  510m. 
The  tubing  is  opened  by  cutting  a  section  out  and  then  winding  the  471  turns 
of  wire  in  the  interior.  After  winding,  the  access  section  is  replaced  and 
the  antenna  sealed.  A  layer  of  7.  62  cm.  wide,  .  013  cm.  thick  copper  tape 
is  wound  and  soldered  around  the  antenna  leaving  a  l/4n  air  gap  at  one 
point.  This  copper  tape  provides  electrostatic  shielding  for  some  of  the 
undesirable  interference.  However,  it  is  not  effective  at  the  lower  ELF 
frequencies  since  a  skin  depth  in  copper  is  greater  than  the  shield  thickness. 
The  shielding  is  effective  at  higher  frequencies  and  is  intended  to  reduce 
strong  E  field  components  from  transmitters  at  VLF  and  above  should  it  be 
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necessary  to  locate  a  measurement  site  near  transmitters.  Such  inter¬ 
ference  can  reduce  the  dynamic  range  of  the  equipment  by  overloading  the 
pre-amplifier. 

Since  the  number  of  turns  are  known  exactly  it  remains  to  establish 

the  loop  area,  from  a  measurement  of  the  mean  loop  diameter.  For  these 

antennas  the  measured  mean  diameter  is  1.  568m  resulting  in  a  turns-area- 

2  . 

product  of  909.  2m  .  Using  the  parameters  developed  here  the  relationship 
between  antenna  voltage  and  field  strength  can  be  rewritten: 

v  =7.  178  x  10'3  |h  If  a-2 

OC  1  0  1 

The  measured  equivalent  circuit  of  the  loop  antenna  is  a  resistor  of 
19.  50  in  series  with  an  inductor  of  0.  8  henries.  The  resistance  is  almost 
entirely  due  to  the  copper  loss  within  the  antenna  since  radiation  resistance 
is  very  small  and  can  be  neglected. 

As  a  check  on  the  performance  of  the  antenna  its  performance  was 
evaluated  using  methods  described  by  NBS[  12]  .  Another  loop  antenna  con¬ 
sisting  of  one  turn  and  a  diameter  of  1.  000  meters  was  constructed  and  used 
in  the  calibration  procedure.  Agreement  between  the  calculated  and  mea¬ 
sured  performance  was  within  the  experimental  accuracy  of  the  measure¬ 
ment  apparatus. 

The  choice  of  a  particular  antenna  configuration  also  involves  a  design 
decision  which  takes  into  account  the  desired  measurements,  ability  to 
match  the  antenna  to  the  pre-amplifier  front  end,  and  the  noise  characteris¬ 
tics  of  the  pre-amplifier.  Originally  the  selection  of  the  1.  568  meter 
diameter  air  core  antenna  was  based  on  a  need  to  make  atmospheric  noise 
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measurements.  In  order  to  maintain  a  20  dB  atmospheric  noise  to  system 
noise  ratio  (1  Hz  BW)  at  20  Hz  the  present  antenna  was  specified  to  realize 
this  ratio  at  that  frequency.  Pre-amplifier  noise  at  these  frequencies  is 
predominantly  l/f  type  and  improved  performance  will  therefore  be  certain  at 
the  higher  frequencies  used  in  the  propagation  tests.  The  point  to  be  made 
here  is  that  the  present  antenna  is  over  specified  for  propagation  measure¬ 
ments. 

A  source  of  noise  is  loop  motion  in  the  earth*s  magnetic  field. 

This  low  frequency  component  associated  with  mechanical  resonances 
in  the  loop  and  mount  can  cause  severe  amplifier  overloading  problems  if 
left  unchecked.  In  most  cases,  the  antenna  has  been  mounted  in  a  hut  to 
isolate  it  from  wind  motion,  and  on  isolation  pads  to  decouple  it  from  earth 
motion  due  to  nearby  vehicular  traffic*  On  some  occasions  the  loops  have 
been  buried  or  partially  buried  in  the  ground  to  avoid  wind  motion. 

A  discussion  of  antennas  is  not  complete  without  mentioning  E  field 

measurements  which  could  be  made  with  a  whip  or  other  type  antenna. 

During  earlier  measurements  of  atmospheric  noise  whip  antennas  were 

employed.  A  great  deal  of  difficulty  was  experienced  in  calibration  of  these 

antennas.  Basically  the  problem  can  be  traced  to  the  high  impedance  levels, 

which  are  present  when  employing  practical  antenna  lengths.  A  source 
7  8 

impedance  of  10  -  10  ohms  is  possible.  At  these  high  impedance  levels 

an  accumulation  of  moisture  from  antenna  to  ground  through  guy  cables, 
base  insulators,  or  the  amplifier  input  terminals  can  seriously  alter  the 
system  calibration.  It  was  found,  for  example,  that  afternoon  rain  in 
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Florida  and  salt  spray  in  Malta  rendered  the  whip  antenna  virtually  useless 
for  dependable  day  in  and  day  out  operation  in  these  areas.  Clouds  at  high 
site  elevations  were  also  a  source  of  error.  Good  results  have  been 
reported  with  electric  field  sensors  from  desert  areas.  In  view  of  the  fact 
that  propagation  measurements  were  planned  world  wide  at  site  locations 
other  than  those  dictated  by  weather  it  was  felt  that  the  whip  antenna  would 
not  produce  the  required  reliability. 

Calibration  Signal  Injection  Methods 

As  shown  in  Fig.  4  the  loop  antenna  output  is  applied  to  a  step-up 
transformer,  pr e- amplifier  and  other  filters  and  amplifiers  in  the  record- 
ing  system.  The  following  is  a  description  of  means  by  which  the  overall 
system  gain  is  calibrated  and  monitored. 

A  calibration  signal  is  continuously  injected  in  series  with  the  antenna 
at  a  known  level  and  frequency.  The  frequency  is  displaced  1  Hz  from  the 
received  signal  frequency.  In  view  of  the  long  recording  periods  (12  hrs.  ) 
which  were  planned  it  was  necessary  to  have  available  a  continuous  record 
of  the  system  performance.  A  frequency  separation  of  1  Hz  was  chosen  as 
a  reasonable  separation  which  would  not  interfere  with  the  received  signal 
analysis  and  yet  would  reveal  possible  changes  in  system  gain  near  the 
frequency  of  interest. 

Referring  again  to  Fig.  4  a  calibration  signal  1  Hz  removed  from  the 
received  signal  is  generated  in  the  synthesizer  and  applied  through  a  preci¬ 
sion  pad  terminated  in  a  0.  500Q  resistor  in  series  with  the  antenna.  A  more 
detailed  schematic  of  the  attenuator  is  shown  in  (Fig.  A-l).  A  measurement 
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Fig.  A-  1 . 


Receiver  front  end 


schematic. 
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of  the  calibration  voltage  at  the  input  of  the  attenuator  pad  is  made  at 
"Cal,  Mon.  M  and  knowing  the  attenuation  ratio  the  voltage  in  series  with  the 
antenna  may  be  calculated.  The  injected  calibration  signal  amplitude  is 
chosen  to  be  sufficiently  large  so  that  it  can  be  accurately  measured  at  the 
output  of  the  receiver.  A  ratio  of  the  detected  calibration  signal  on  playback 
to  injected  signal  at  the  antenna  gives  the  gain  of  the  system  at  the  calibra¬ 
tion  frequency. 

It  will  be  noticed  that  the  attenuator  pad  is  of  the  "balanced  H"  type. 

It  is  made  up  of  .  05^  and  .  01^  resistors  and  has  a  calculated  attenuation  of 

-  7 

4.  4959  x  10  .  The  attenuation  of  this  pad  has  been  measured  in  the  labora¬ 

tory  and  agrees  within  experimental  error  with  the  calculated  value.  Since 
the  accuracy  of  the  test  equipment  used  in  this  measurement  is  less  than 
the  combined  accuracy  of  the  resistors  the  calculated  value  of  attenuation 
has  been  utilized  throughout.  Although  the  attenuator  pad  itself  is  balanced 
it  is  employed  in  an  unbalanced  circuit  since  one  side  of  the  pad  is  grounded 
along  with  pins  1,  8  and  9  of  the  input  transformer.  A  short  discussion  of 
the  reasons  behind  this  arrangement  follows. 

Initial  attempts  to  employ  a  balanced  input  revealed  a  common  mode 
voltage  (with  respect  to  ground)  on  the  balanced  pad.  This  signal  was  intro¬ 
duced  by  the  capacitive  coupling  between  windings  on  the  calibration  input 
transformers.  This  unwanted  signal  was  capacitively  coupled  across  the 
input  transformer  windings  to  the  unbalanced  input  of  the  40  dB  Ithaco  pre¬ 
amplifier.  As  such^it  added  vectorially  with  the  desired  injected  signal 
which  was  coupled  through  the  transformer,  and  gave  an  erroneous  measure 
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of  the  system  gain.  (Note  that  the  relatively  high  impedance  levels  in  the 
circuits  make  it  relatively  easy  for  11  stray"  capacitances  to  inject  unwanted 
signals.  )  Various  configurations  of  balanced  and  unbalanced  inputs  were 
tried  with  the  method  shown  giving  the  best  results.  Measurements  show 
that  the  common-mode  signal  is  30  -  40  dB  below  the  desired  signal  with 
this  arrangement. 

A  calibration  correction  factor  is  computed  which  compensated  for 
differential  system  gain  at  the  calibration  and  signal  frequencies.  Differ¬ 
ences  do  exist  due  to  ripple  in  band  pass  filters,  mis-tuning  of  filters,  and 
drift  of  component  values  with  time  and  temperature.  In  order  to  mitigate 
this  problem  the  system  gain  is  measured  with  a  shorted  antenna  at  the  two 
frequencies  of  interest  before  and  after  a  recording  session.  This  informa¬ 
tion  is  entered  into  the  log  and  later  used  as  a  correction  factor  in  the 
signal  estimate.  It  has  been  found  from  experience  that  the  gain  at  the  two 
frequencies  is  nearly  the  same  and  does  not  change  materially  over  a 
measurement  session  (3-4  weeks).  Of  course,  it  is  possible  to  experience 
a  system  gain  change  over  several  hours  but  this  will  be  noticeable  in  the 
calibration  level  on  playback  and  is  correctable. 

The  input  to  the  attenuator  pad  is  measured  with  a  Ballantine  model 
3Z3  True  RMS  Voltmeter.  It  has  a  rated  accuracy  of  Z 4>  at  all  points  on 
scale.  Before  each  measurement  period  all  meters  used  in  the  experiment 
are  compared  with  a  laboratory  secondary  standard  at  the  level  where  they 
will  be  used  (110  -  1Z  5  mv).  From  this  it  is  felt  that  the  meters  are  accurate 
to  within  0.  1  $  absolutely  and  have  negligible  error  differentially.  The 
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largest  source  of  error  in  the  field  is  probably  due  to  meter  reading. 

Front  End  Considerations 

Perhaps  the  most  important  design  problems  are  encountered  in  the 
front  end.  A  good  design  will  integrate  antenna  characteristics  with  the 
pre-amplifier  to  realize  a  low  noise  front  end.  Originally  the  front  end 
design  was  based  on  a  need  to  measure  wide-band  atmospheric  noise  with 
high  fidelity.  As  an  interest  in  both  propagation  and  atmospheric  measure¬ 
ments  developed  the  front  end  design  was  reevaluated  and  found  to  satisfy 
both  requirements. 

A  great  deal  of  early  effort  was  expended  in  identifying  a  low-noise 
amplifier  suitable  for  use  at  the  frequencies  and  impedance  levels  which 
were  present.  Source  impedances  would  be  of  the  order  of  1  megohm  or 
greater  and  the  range  of  frequencies  extended  from  20  -  300  Hz.  Under 
these  conditions  the  Field  Effect  Transistor  (FET)  is  well  suited  for  ampli¬ 
fier  use.  Figure  A-2  shows  the  equivalent  noise  sources,  both  voltage  and 
current  generators,  which  are  useful  in  characterizing  the  noise  perform¬ 
ance  of  the  amplifier.  A  recent  tutorial  article  [  13]  describes  the  uses  of 
this  equivalent  circuit,  and  follows  the  design  procedure  employed.  A  major 
result  of  this  analysis  is  the  relationship  for  the  optimum  source  impedance 

for  low  noise  operation.  Optimum  source  impedance  is  given  by: 

e 

Z  =  . 

opt  1 
r  n 

Since  e^  and  i  are  both  functions  of  frequency  the  optimum  source  impedance 
varies  across  the  band  of  interest  and  is  difficult  to  match  over  a  wide  band. 
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Fig.  A-2.  Equivalent  circuit. 
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Several  commercial  FET  amplifiers  were  investigated  in  the  prelimi¬ 
nary  design.  In  addition  a  laboratory  designed  amplifier  was  evaluated.  As 
a  result  of  these  investigations  it  was  concluded  that  all  reasonably  compe¬ 
tent  amplifier  designs  resulted  in  performance  characteristics  which  were 
very  similar.  It  was  also  concluded  that  single-ended  input  amplifiers 
resulted  in  lowest  noise  operation.  Parenthetically  it  might  be  noted  that  a 
selection  process  is  employed  for  the  crucial  input  transistor  and  each 
supplier  had  his  own  proprietary  circuit  and  method.  The  amplifier  that 
was  selected  as  a  result  of  this  investigation  was  the  Ithaco  Model  144.  This 
amplifier  meets  all  of  the  requirements  for  the  system  and  has  the  most 
detailed  specifications  of  all  of  the  units  examined. 

Referring  to  the  Ithaco  specification  sheet  the  equivalent  noise  current 
in  a  1  Hz  bandwidth  is  typically  .  01  picoampere  in  the  10  to  100  Hz  range. 

-  8 

Shorted  source  equivalent  noise  voltage  in  the  same  bandwidth  is  Z.  5  x  10 
volts  at  50  Hz.  Utilizing  these  two  pieces  of  information  the  optimum  source 
resistance  at  50  Hz  is  calculated  as: 

%  =  2.5x_10-8  .  ,  5  x  106n 

opt  »n  10-4 

Note  that  due  to  the  l/f  behavior  of  the  equivalent  voltage  noise  source,  the 
optimum  source  impedance  will  decrease  with  frequency. 

In  order  to  obtain  a  match  between  the  antenna  and  the  amplifier  to 
effect  low  noise  performance  a  50:1  step-up  transformer  is  employed  be¬ 
tween  the  antenna  and  the  pre-amplifier  .  At  100  Hz  the  antenna  impedance 
is  given  by: 
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=  19.  5  +  j502Q  N  j502n 


Zant 

Neglecting  transformer  primary  winding  resistance  and  secondary  winding 
capacitance  the  antenna  equivalent  impedance  reflected  to  the  secondary  is: 

Z  «  (50)2  (j502)  =  jl.  25  X  106n. 

source  '  '  9  J 

A  2.  2  x  106n  resistor  is  placed  across  the  transformer  winding  to  flatten 
the  frequency  response  of  the  front  end.  Including  this  resistance  the  mag¬ 
nitude  of  source  impedance  is 

Iz  I  =  1. 1  x  io6n. 

1  sourc  e 

It  is  observed  that  the  optimum  source  impedance  is  not  achieved  at 
this  frequency;  however,  the  degradation  due  to  this  will  be  small.  Due  to 
the  inductive  nature  of  the  antenna  its  contribution  to  the  source  impedance 
will  increase  with  frequency  whereas  the  optimum  impedance  should  fall  off. 
In  the  frequency  band  40  -  300  Hz  front  end  noise  is  at  least  30  dB  below 
median  atmospheric  noise  with  this  arrangement. 

In  the  equipments  employed  for  atmospheric  noise  measurements  the 
front  end  was  located  adjacent  to  the  antenna  and  power  was  supplied  from  a 
separate  battery  supply.  After  some  amplification  the  received  signal  was 
transferred  via  a  balanced  output  line  driver  over  a  twisted  pair  to  the  rest 
of  the  recording  equipment.  This  arrangement  meant  that  the  front  end 
would  be  exposed  to  the  extremes  of  temperatures  and  humidity.  In  fact, 
oscillations  did  develop  with  some  equipments  in  the  presence  of  high  humi¬ 
dity.  Excessive  moisture  formed  feedback  paths  to  the  pre-amplifier  input 
via  the  printed  circuit  boards  that  were  employed.  The  present  equipments 
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utilize  a  front  end  located  in  a  rack  with  the  remainder  of  the  equipments. 
Coupling  to  the  antenna  is  accomplished  via  No.  22  shielded  twisted  pair. 
Twisting  avoids  incidental  magnetic  pick-up  and  the  shielding  reduces  the 
effects  of  electric  fields.  The  antenna  is  located  approximately  200'  from 
the  recording  equipment  and  300f  from  the  DC  motor  generator  to  avoid 
magnetically  generated  interference. 

Several  other  precautions  have  been  taken  to  avoid  interference  from 
other  parts  of  the  system  and  these  will  be  described  briefly  here.  The 
front  end  including  transformer,  amplifier  and  low-pass  filter  are  enclosed 
in  a  50  mil  high  nickel  steel  alloy  box  to  avoid  magnetic  coupling  from  other 
components.  Such  shielding  is  effective  against  local  power  line  interfer¬ 
ence  and  possible  coupling  with  the  locally  generated  reference.  This 
shielding  is  in  addition  to  that  supplied  by  the  transformer.  Another  import¬ 
ant  consideration  is  the  ground  system.  All  grounds  are  made  with  heavy 
copper  braid  and  are  common  at  one  point.  Ground  return  currents  which 
might  induce  unwanted  voltages  at  the  amplifier  input  are  avoided.  Although 
the  Ithaco  amplifier  exhibits  55  dB  of  power  supply  isolation  great  care  is 
taken  to  de-couple  the  power  supply  voltage  which  feeds  the  amplifier.  If 
suitable  de-coupling  is  not  employed,  the  amplifier  noise  performance  may 
be  degraded,  or  more  seriously,  components  of  the  locally  generated 
reference  may  be  coupled  into  the  front  end.  Finally,  the  front  end  and  all 
other  components  of  the  recording  system  are  subject  to  environmental 
control. 
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Notch  Filters 


Notch  or  band  reject  filters  are  included  in  the  signal  path  to  reduce 
unwanted  pick-up  from  power  lines.  Several  card  slot  positions  are  alloca¬ 
ted  for  this  purpose  in  order  that  the  most  effective  combination  of  filters 
may  be  synthesized  to  combat  the  problems  at  the  particular  site.  A  normal 
complement  of  filters  in  CONUS  would  include  on  each  of  60,  180  and  300  Hz 
rejects.  Overseas  operation  at  sites  like  Norway  or  Malta  would  utilize 
50  Hz  and  higher  order  harmonic  rejection.  Experience  has  shown  that  the 
fundamental  and  odd  order  harmonics  cause  the  severest  problem.  In  fact, 
the  Malta  site  has  had  the  dubious  distinction  of  having  a  larger  component 
at  150  Hz  than  50  Hz. 

Another  problem  that  is  sometimes  encountered  is  that  of  frequency 
wander  in  the  local  power  system.  Frequency  variations  of  ±0.  5  Hz  have 
been  observed  in  Malta.  It  has  been  profitable  in  this  situation  to  stagger 
tune  three  notch  filters  to  49.  5,  50  and  50.  5  Hz  to  establish  band  rejection 
from  49.  5  to  50.  5  Hz. 

Individual  notch  filters  have  been  designed  using  methods  described 
in  a  Lincoln  Laboratory  Technical  Note  [  14  ]  •  Typically  the  notch  filters 
exhibit  40  dB  of  rejection  at  their  center  frequency.  Quality  factor  ("Qn)  is 
adjusted  to  be  in  the  range  from  8  to  10.  All  components  are  measured 
within  0.  1  $  before  assembly,  and  final  adjustments  are  made  on  each  filter 
using  a  precise  frequency  source.  Where  several  filters  are  cascaded  to 
obtain  a  wider  rejection  range,  the  individual  filters  are  adjusted  in  fre¬ 
quency  and  MQM  to  obtain  the  desired  response.  Insertion  loss,  out  of  band, 
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is  typically  6  dB. 
Band  Pass  Filters 


Band  pass  filters,  sometimes  referred  to  as  narrow  band  filters,  are 
employed  in  the  signal  channel  to  establish  a  narrow  band  channel  centered 
on  the  receiver  frequency.  Narrowing  the  bandwidth  prior  to  recording 
eliminates  problems  associated  with  the  large  dynamic  range  of  wide  band 
atmospheric  noise  as  opposed  to  the  limited  dynamic  range  of  the  tape 
recorder  (40  dB).  All  signal  estimates  are  made  from  the  narrow-band 
channel  using  linear  processing.  Individual  filters  are  plugged  in  as  opera¬ 
tions  demand,  each  filter  covers  three  or  more  possible  transmitter  fre¬ 
quencies. 

The  narrow  band  filters  consist  of  two  single  pole  stagger  tuned  band 
pass  filters  mounted  on  a  single  card.  Typical  bandwidths  in  the  40  -  50  Hz 
range  are  3  Hz,  (and  4  Hz  in  the  70  -  80  Hz  range).  Both  11  Qn  and  center 
frequency  of  the  individual  filters  are  adjusted  to  obtain  desired  response. 
Insertion  gain  varies  from  24  -  28  dB  depending  on  center  frequency. 

Design  procedures  are  the  same  as  those  employed  for  the  notch  filters 
[  14  ]. 

Required  Frequency  Stability 

In  view  of  the  long  integration  periods  (1-2  hours)  used  to  estimate 
the  received  signal  it  is  necessary  to  examine  the  required  frequency  accu¬ 
racy  to  avoid  degradation  of  the  estimate.  The  model  shown  in  Fig.  A- 3  is 
representative  of  the  measurement  process  and  will  be  used  to  derive  the 
required  frequency  stability. 
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Figure  A-3 

The  imperfect  local  reference  can  be  written  as: 

cos  it)  t  cos  U)  0 1  -  sin  up  t  sin  (13  et 
s  s  s  s 

The  output  of  the  multiplier  is  given  by: 

2 

A  (cos  uu  t  cos  UP  et  -  sin  up  t  cos  up  t  sin  up  et) 
s  s  s  s  s 


=  {  [  cos  Zup  t  +  1]  cos  U)  Gt  -  [  sin  2tu  t]  sin  UP  et}  . 

C*  s  s  s  s 

After  removal  of  double  frequency  terms  with  the  low  pass  filter  the  inte¬ 
grator  output  is: 

A  T  A  sina,seT 

eo  *  T  J  „  cos  V*  dt  =  "T  i  t 

U  s 

where  T  is  the  integration  period. 

Now  expanding  sin  up  GT  into  a  series  and  keeping  the  first  two  terms  yields, 


e  =  ■ 


o  UP  0T 
s 


(ai  GT y 

[i^s^T  -  — - ] 


when  the  second  term  is  much  smaller  than  the  first  then  e  =  A  as  desired. 

o 

The  frequency  stability  for  a  l4>  error  (.  1  dB)  is  then  computed  from 

,2 


(ouseT)‘ 
TT 


=  .  01 


solving  for  e 


94 


Using  worst  case  parameters, 

_3 

T  =  2  hrs  =  7.  2  x  10  seconds 


fs  =  80  Hz 


.  245 


6.  75  x  10 


-8 


G 


- T" 

(zn  X  80  X  7.  2  X  10  ) 


That  is,  a  frequency  difference  between  transmitter  and  receiver  of 

g 

6.  75  parts  in  10  will  cause  a  0.  1  dB  degradation  in  signal  estimation  over 
a  Z  hour  integration  period. 

A  short  discussion  of  possible  frequency  errors  is  now  in  order. 

Prior  to  deployment  in  the  field  the  receiver  frequency  standards  are  com¬ 
pared  to  a  Lincoln  Laboratory  Cesium  Beam  Standard  which  is  itself  tied 
to  NBS.  The  transmitter's  frequency  standard  is  also  continuously  com¬ 
pared  to  NBS  (WWVL)  and  maintained  to  within  1  part  in  10^.  Initially,  the 
receiver  frequency  standards  are  set  to  within  2  parts  in  10^  initially,  and 
typically  exhibit  a  worst  case  monotonical  change  in  frequency  of  2  parts 
in  1010  per  day.  Approximately  320  days  would  have  to  elapse  before  a 

g 

discrepancy  of  6.  75  parts  in  10  was  reached.  From  this  it  can  be  seen 
that  errors  due  to  frequency  drift  have  been  negligible. 

Frequency  Generation 

Two  precise  frequencies  are  generated  and  recorded  at  each  site. 
These  two  frequencies  have  been  entitled  the  "Calibration"  and  "Signal" 
reference  frequencies.  The  calibration  signal  is  1  Hz  removed  from  the 


95 


received  signal  and  is  recorded  as  a  sine  wave  at  high  level  on  channel  7. 

An  attenuated  version  of  this  signal  is  inserted  in  series  with  the  antenna 
and  is  employed  on  playback  to  ascertain  system  gain.  The  Signal  reference 
is  a  square  wave  at  the  received  signal  frequency,  and  is  recorded  at  a  low 
to  intermediate  level  on  channel  6.  Both  a  frequency  standard  and  a  synthe¬ 
sizer  are  required  to  generate  the  required  frequencies. 

The  frequency  standard  is  a  Frequency  Electronics  Model  1000-A  with 
a  frequency  stability  of  1  -  Z  parts  in  10^  per  day.  A  battery  pack  is  pro¬ 
vided  to  allow  8  -  Z4  hours  operation  in  the  event  of  primary  power  loss. 
Power  for  the  frequency  standard  when  operating  in  the  system  is  supplied 
at  4-30  VDC. 

Prior  to  conducting  measurements  at  field  sites  all  frequency  stand¬ 
ards  are  compared  with  a  Cesium  Beam  Standard  located  at  the  laboratory. 
Over  a  several  day  period  all  standards  are  brought  into  agreement 
(~  Z  parts  in  10^)  with  the  Cesium  Beam  Standard  and  their  behavior  moni¬ 
tored.  All  standards  are  under  power  continuously  and  travel  with  site 
personnel  to  and  from  the  sites.  If  a  transit  time  of  more  than  8  hours  is 
required  an  additional  battery  pack  is  f,floatedM  across  the  internal  supply 
to  assure  continuous  operation.  Whenever  possible  the  frequency  accuracy 
of  the  standards  is  compared  against  other  systems  in  the  field.  For  this 
purpose  comparisons  with  WWVL  (60  kHz)  using  the  HP  model  117A  VLF 
Comparator  and  comparisons  with  LORAN  (100  kHz)  using  the  Beukers 
Model  11Z  Frequency  Reference  have  been  used.  Upon  completion  of  the 
experiment  the  standards  are  returned  under  power  to  the  laboratory  and 
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their  frequency  again  compared  to  the  Cesium  Standard.  A  back-up  standard 
is  available  at  each  site  in  the  form  of  the  1  mHz  crystal  oscillator  in  the 
Time  Code  Generator  (TCG).  Standard  operating  procedure  is  to  compare 
the  FEL  standard  with  the  TCG  on  site  and  adjust  as  closely  as  possible  the 
TCG  oscillator  so  that  it  agrees  with  the  standard.  If  the  standard  should 
fail  it  is  then  possible  to  continue  with  the  back-up  oscillator  for  a  limited 
time. 

The  synthesizer  accepts  the  1.  0  mHz  output  from  the  Frequency 
Standard  and  generates  the  two  required  frequencies  using  two  independently 
controlled  sub- synthe sizer s  with  identical  logic.  Output  of  the  synthesizer 
consists  of  two  independent  frequencies  selectable  in  1.  0  Hz  steps  from  1 
to  1Z7  Hz.  Ordinarily  frequencies  between  40  to  50  Hz  and  70  to  80  Hz  are 
employed  since  this  corresponds  to  the  present  transmitter  capability.  The 
philosophy  behind  the  synthesizer’s  operation  can  be  explained  briefly  as: 

1.  A  high  repetition  rate  pulse  train  (average  frequency  960  kHz) 
is  fed  into  a  six  stage  counter.  A  non-overlapping  transition  from  each 
stage  plus  the  input  is  used  to  generate  a  narrow  pulse  thus  establishing 
seven  pulse  trains  related  in  frequency  by  powers  of  2. 

2.  Combinations  of  these  pulse  trains  are  selected  by  front  panel 
switches  and  combined  in  an  f,ORlf  gate. 

3.  Additional  countdown  circuits  cause  the  output  to  appear  in  the 
frequency  range  of  interest. 

4.  The  calibration  frequency  reference  is  applied  to  a  low  pass 
filter  and  appears  as  a  sine  wave.  The  signal  frequency  reference  is 
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employed  as  a  square  wave. 

Notice  that  this  method  makes  use  of  the  average  repetition  rate  of  the 
higher  frequency  pulse  train,  the  consequence  of  which  is  a  waveform  with 
jitter.  After  further  division  by  a  large  number  this  jitter  is  reduced  pro¬ 
portionally  and  is  typically  10  psec  at  the  output. 

Circuitry  to  detect  phase  jumps  in  the  signal  reference  frequency  has 
been  incorporated  in  the  synthesizer.  The  purpose  of  this  circuitry  is  to 
detect  possible  phase  jumps  due  to  faulty  operation  either  of  the  logic  or 
external  causes.  If  such  a  jump  is  detected  a  light  is  activated  which  may 
only  be  extinguished  manually.  Since  the  equipment  is  required  to  operate 
for  1Z  hour  intervals  in  an  unattended  condition  it  was  felt  that  brief  interrup¬ 
tions  of  the  reference  could  occur  and  go  unnoticed  unless  some  precaution 
was  taken.  External  influences  such  as  nearly  mobile  radios,  and  starting 
of  vehicles  has  caused  phase  jumps  in  the  past.  Nearby  thunderstorms 
such  as  those  experienced  in  Florida  could  also  be  expected  to  trigger  the 
logic  occasionally. 

Tape  Recorder 

Recording  at  the  field  sites  is  accomplished  with  a  Precision  Instru¬ 
ment  Model  P15100  l/zM  tape  recorder.  Tape  speed  is  15/16  IPS  and 
10  -  l/Z"  reels  with  1.  0  mil  mylar  tape  is  utilized.  This  combination 
supplies  a  1Z  hour  record  time.  Seven  high  quality  FM  data  channels  with 
approximately  31Z  Hz  bandwidth  are  used  as  the  primary  data  channels. 

Two  back-up  lower  quality  edge  track  channels  are  also  available,  but  were 
not  employed.  Power  consumption  is  approximately  18  watts  including 
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playback  monitor.  Tape  speed,  wow  and  flutter  are  specified  as  ±  1$  or 
less. 

Since  any  analog  recording  system  exhibits  a  tape  speed  variation 
with  time  it  is  necessary  to  record  reference  signals  at  the  time  of  record¬ 
ing  if  long  term  integration  is  required  on  playback.  In  order  to  support 
the  measurement  requirements  of  the  system,  a  square  wave  at  the  same 
frequency  as  the  received  signal  is  recorded  on  channel  6,  and  a  high  level 
calibration  signal  is  recorded  on  channel  7.  These  signals  are  used  as 
references  on  playback  and  serve  to  eliminate  tape  speed  variations  which 
occur  on  record  or  playback. 

The  tape  recorders  have  demonstrated  a  degradation  in  performance 
over  a  three  week  measurement  period.  This  degradation  is  evidenced  by 
variations  in  tape  speed  of  -fc  5$  over  a  12  hour  interval  and  increased 
flutter  particularly  towards  the  end  of  a  12  hour  tape.  A  large  portion  of 
these  troubles  has  been  traced  to  a  carbon  disc  slip  clutch  in  the  feed  reel 
as  well  as  improper  adjustment  of  pinch  rollers.  Since  the  original  set  of 
measurements  the  tape  recorders  have  been  returned  to  the  factory,  modi¬ 
fications  made,  and  improved  line-up  procedures  developed.  It  is  extre¬ 
mely  important  that  the  manual  accompanying  the  tape  recorder  be  used  as 
a  reference  for  all  adjustments  and  maintenance,  and  that  a  complete 
checkup  be  conducted  prior  to  field  use. 

Prime  Power 

Prime  power  is  supplied  by  an  Onan  AK  series  DC  motor  generator 
set  wdth  a  capacity  of  750  watts.  The  motor -generator  set  is  operated 
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continuously  except  for  oil  changes,  and  supplies  a  charging  current  for  a 
set  of  high  capacity  batteries  (70  -  100  amp-hour)  which  in  turn  supply  the 
required  current  at  28  volts  to  the  recording  system.  A  modification  to  the 
field  winding  of  the  generator  allows  the  generator  output  to  be  balanced 
against  the  system  demand  with  a  net  trickle  current  (.  5a)  available  for  the 
batteries.  Battery  capacity  is  sufficient  to  allow  proper  operation  of  the 
site  for  a  period  of  12  hours  should  the  motor-generator  set  fail.  A  twelve 
gallon  tank  capacity  is  provided  and  the  tank  is  replenished  daily  with 
regular  grade  gasoline. 

Care  is  taken  when  setting  up  a  site  to  locate  the  motor  generator  set 
as  far  away  as  possible  from  the  antenna  and  recording  gear.  This  arrange¬ 
ment  tends  to  eliminate  both  electrical  and  mechanical  coupling  to  the 
antenna.  Electrical  interference  due  to  the  motor-generator  set  has  never 
been  noticeable  unless  the  antenna  and  generator  were  within  50  feet  of 
each  other.  Mechanically  coupled  interference  due  to  induced  motion  of  the 
antenna  in  the  earth's  magnetic  field  has  been  observed  at  distances  up  to 
250  feet  in  the  situation  where  the  antenna  was  buried  and  the  motor -genera¬ 
tor  was  resting  on  the  ground.  A  simple  solution  to  this  problem  was  to 
place  the  motor  generator  on  several  old  tire  carcasses  to  insure  mecha¬ 
nical  de-coupling  from  the  ground. 

If  the  maintenance  suggestions  in  the  Onan  Manual  are  followed  the 
generator  should  give  satisfactory  service.  The  following  is  a  short  list 
of  some  of  the  more  important  precautions: 

1.  Use  regular  grade  gasoline 


J00 


2.  Change  the  oil  weekly  using  MS  or  MS/DS  service,  preferably 
Shell 

3.  Remove  the  panels  covering  M.  G.  set  in  a  hot  climate,  and  at 
altitudes  exceeding  4000  feet 

4.  Re-adjust  the  carburetor  towards  the  lean  side  when  operating 
at  high  altitudes 

5.  Dismantle  the  engine  and  clean  the  carbon  after  500  hours. 
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APPENDIX  B 


Noise  Intensity  Measurements 

Accurate  knowledge  of  the  noise  intensity  (at  the  signal  frequency)  is 

necessary  for  the  following  two  reasons. 

1.  The  estimate  of  the  signal  power  at  the  output  of  the  receiver's 

^2  2  No 

summer  is  a  biased  estimator  (i.e.,  E[A  ]  =  A  +2  (—^r  )) .  The  bias  which 
is  determined  by  the  noise  intensity  consequently  must  be  subtracted  from 

A 

the  sample  value  A  to  give  the  unbiased  estimate  of  A  . 

A 

A  2.  The  confidence  limits  on  each  estimate  A  depend  on  the  ratio 

2  N  A 

A^  1  o  .  2 

(where  N^  =  -7^-).  Accordingly  a  small  confidence  interval  for  A 

requires  an  even  smaller  confidence  interval  for  N^. 

For  purposes  of  this  experiment  we  have  set  an  acceptable  confidence 

interval  of  ±  .  2  dB  on  Nj.  It  should  be  noted  that  the  best  way  to  reduce  the 

confidence  interval  in  estimating  a  signal  level  is  to  integrate  over  the  entire 

A 

data  record  and  come  up  with  a  single  value  A  .  On  the  other  hand,  to 

reduce  the  confidence  interval  in  noise  intensity  estimation  the  scheme  is  to 

get  as  many  independent  samples  in  the  data  record  as  possible.  We  will 

defer  until  the  end  of  this  section,  the  method  by  which  we  select  the  number 

of  samples  needed  to  provide  a  specified  confidence  interval  about  N^.. 

Figure  B- 1  shows  the  primary  method  by  which  the  noise  intensity  was 

computed.  The  output  of  a  single  correlator  (PAR)  is  sampled  at  a  reasonably 

high  rate  (^gamp^e  ~  20/sec)  and  recorded  on  magnetic  tape.  The  output 

random  variable  is  denoted  as  L  (x  denotes  either  the  cosine  or  sine 

x 

—  2 

reference  channel).  The  sample  mean  m  and  sample  variance  q-  are 

X  X 
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defined  by 


(b-  1 ) 


_  1  M 

m  =  Z  (L  ). 

x  M  .  ,  x  i 

i=  1 


2  _  J_  ^  r  /  T  \  - 

°x  M2,  x  i  mx^ 
1=  1 


(b-2 ) 


It  can  be  shown 


KA  A  _  KAiA2 

E[m  ]  =  ( - -)  A  cos  (p  ,  E[m  ]  =  ( — — - -)  A  sin  (n 

c  Jz  Jz 


and 


-  ,  KA.A,  2  N  KA.A  2  N  .  N'  , 

r  2  _  ,M- 1  12,  o  ,  _  12,  o^  A  o,l 

Elax  1  -  (— >  —  -  (  T - ~ 


(b-3) 


for  large  M  where  N  =  N  n  in  Fig.  b-1;  Nf  ,  includes  the  constant 
b  o  o,  1  o,  1 

(KAjA^/v^Z)^  and  T  =  2(32)RC.  (RC  being  the  time  constant  of  the  LPF.  ) 
The  above  is  true  provided  the  random  variables  (L»x)^  are  identically 
distributed  and  statistically  independent. 


The  important  point  is  that  the  sample  variance  (o^  )  is  computed 
from  the  M  independent  samples  taken  at  the  output  of  one  of  the  correlators. 
The  routines  for  sample  mean  and  variance  have  been  implemented  on  a 
Varian  620i  computer.  Since  the  received  signal  levels  are  expressed  in 

terms  of  an  H  field  at  the  input  to  the  antenna,  we  need  to  relate  an 

2  ~ 
equivalent  input  noise  intensity  to  the  Varian  computed  variance  crx  . 

First  the  noise  intensity  at  the  antenna  output  [  (volts)  /Hz]  is  related  to  an 

equivalent  magnetic  field  noise  intensity  [  (amp-m)  /hz]  though  equation  a-2 

N 

^  ;=»  m  r 

(b-4) 


H 


o,  ant 


N  (.007178f)2 


The  system  gain  at  the  signal  frequency  from  the  antenna  output  to  the  tape 
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recorder  output  is  the  product  of  tape  recorder  voltage  gain  G,^  and  the 

voltage  gain  between  the  antenna  output  and  tape  recorder  input 

G  G  V  (f  ) 

Gsys  =  Gt,GS  =  GTGC(G^)  =  GT(G^)(  g°f  ')  }  ’  (b"5) 

The  parameters  in  equation  b-5  are  fully  defined  in  Appendix  E  and  it  is  only 
necessary  to  note  here  that  they  are  all  quantities  measured  during  the 
normal  receiver  calibrations.  Along  with  the  tape  recorder  voltage  gain 
there  is  a  noise  bandwidth  expansion  due  to  playback  speedup.  The  effect  of 
speedup  is  to  stretch  the  noise  spectrum  by  the  speedup  factor  F.  Follow¬ 
ing  the  tape  recorder  we  have  been  using  an  attenuator  (with  voltage  gain 

<  1),  Hence  the  noise  intensity  input  into  the  (PAR)  correlator  is  related 

to  the  antenna  output  noise  intensity  by 

N  .  .  (G  .  P,)2 

N  ,  =  °’ant  -  s?s - l—  .  <b-6> 

o,  1  F 


The  output  of  the  PAR  correlator  (Nq  can  be  shown  to  be  related  to  the 
input  noise  intensity  Nq  ^  by  the  relationship 


N  l  I  2 

No,  4  =  2&T  ’  • 


(b-7) 


Finally,  the  voltage  amplifier  A^  and  the  A/D  voltage  gain  [A^]  relate 


^  and  the  computed  variance 


»x2  ■  <AlV2No,4  • 

Equations  b-4,  b-5,  b-6,  b-7  and  b-8  are  combined  to  yield 


(b-8) 


H 


N  (.007l78f)2  (Gt.Gs.P1)2 


4RC 

K2 


(aia2) 


2  CTx 


(b-9) 


For  the  actual  analysis  the  following  constants  were  used  A^  =  A^=  1, 
RC=  .  0 1  sec,  F=  32  and  K=(10/(50  x  10'3)). 
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H^  is  the  equivalent  antenna  input  noise  intensity  f(amp-m)  /Hz)  of  noise 
picked  up  by  the  single  loop  antenna.  The  results  of  the  noise  measurements 
are  listed  in  tables  VIII  and  IX  in  Section  V. 

The  virtue  of  the  above  noise  analysis  is  that  a  number  of  noise 
estimates  for  each  tape  can  be  computed  several  minutes  after  the  equipment 
has  been  set  up  and  calibrated.  An  alternative  method  of  computed  noise 
intensity  is  used  on  occasion.  The  method  consists  of  digitizing  large 
segments  of  the  wideband  channel  on  playback.  The  resulting  digital  tape  will 
be  in  the  same  format  as  was  used  for  Lincoln's  wideband  noise  study.  One 
of  the  spectral  programs  produces  a  noise  spectrum  from  which  one  can 
pick  off  the  noise  intensity  level  at  45  Hz.  The  primary  disadvantage  of  this 
method  is  that  it  is  time  consuming  and  costly  in  computer  time.  The 
agreement  between  the  noise  intensity  made  by  variance  computation  and  that 
made  by  the  spectrum  estimate  is  extremely  close  (differences  are  seldom 
more  than  a  few  tenths  of  a  dB).  The  advantage  of  the  spectrum  computation 
is  that  it  shows  clearly  sources  of  strong  interference  at  frequencies 
sufficiently  close  to  45  Hz  (7  5  Hz)  which  may  affect  the  signal  estimate.  The 
spectrum  computations  also  clearly  show  the  noise  to  be  white  over  the 
receiver  bandwidths  of  interest. 

As  a  slight  digression  we  will  outline  how  we  select  the  number  of 
samples  needed  to  satisfy  a  specified  confidence  interval  about  H^(dBH). 

Let  us  first  assume  that  the  output  of  the  correlator  (L^).  has  had  its  mean 
subtracted  off.  Denote  this  new  variable  by  Previous  empirical  results 
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derived  from  Lincoln's  wideband  noise  studies  of  ELF  noise  have  indicated 
that  as  the  non-Gaussian  wideband  noise  is  subjected  to  successively 
narrower  band  pass  filters,  the  output  becomes  more  Gaussian.  This  is 
apparent  for  real  time  equivalent  integration  times  of  a  few  seconds  for  most 
sites  where  wideband  noise  has  been  collected.  The  probability  distribution 
function  for  L  is  then  approximated  by 


%<v 


(b-10) 


The  sample  variance  (equation  b-2)  can  be  written 


2  1  ^  ,  2  A  1  2 

°  ~  M  li  ~  M  X 

1=  1 

0^  itself  is  then  a  random  variable  with  a  distribution  which  can  be  shown  to 


be  a  x  distribution  with  M  degrees  of  freedom,  i.  e.  , 

f  i  -  M  ->M/2  ,M.  r  i 

f  2(ct  *  (N 7  ,/T)  Z  r  (2  >  (N1  ,/T)  ^ 

0*  ®  9  O  i  A 


»  r  Mpr _ 

2  f  -  1  *(n£  ,/T) 

X  e 


(b-  11) 


for  o'  ^  0 
=  0  for  aZ  <  0. 

The  mean  and  variance  of  the  random  variable  o^  are  most  easily  found  by 
determining  the  moment  generating  function  for  the  random  variable  . 
This  can  be  shown  to  be 

^  2(t) 


J—W 7 


N 1  , 

2  (-§4) 
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Since  the  random  variables  -F  are  sampled  at  a  sufficiently  slow  rate  so  as 

to  imply  independence  of  samples  one  can  construct  the  moment  generating 

A  M  2 

function  of  the  sum  Z  =  £  l.  from 

i=l  1 

M 

^  (t)  =  0  (t)  =  H  ib  ?(t) 

^  Zl.  i=  1  l 

1  1 


i/j  ^  (t)  = 


l 


N 


(2- 


'  '  M/2 

O.  Is  ' 


2  ( 


N'  , 

o,l, 


-  t} 


M/2 


(b-  12) 


Now  the  first  two  moments  of  the  distribution  of  Z  are  found  from  equation 
b-  12 

?  i  i  9  7  (t)  1 

(b  -  13) 


p,  r  2  1  1  90Z(t) 

E[a  ]  =  -  E[Z]  =  - 


a  t 


N 1  , 

I  =  (  > 

1 1=  0  '  T  1 


and 


var[(TZ]  =  -L  {E[ZZ]  -  (E[Z])Z}  =  |  -  |  )Z  }. 

dr  t_u  dr  l~ 


N'  ,  2 

2(-4lJ-) 


M 

2  No  1 

From  b-13  one  sees  that  ^  is  an  estimator  of  — ^ —  and  from  b-14  the 
variance  of  the  estimate  decreases  inversely  with  the  number  of  samples 


(b-14) 


chosen.  It  is  observed  that  f  )  is  a  function  both  of  M  and  the  true 

No  1  °  2  2 

variance  ( - — — )  •  The  confidence  limits  (7tt  and  are  chosen  such  that 


+  00 


N1 

Oj_L,  a  ^  _  _£ 


J  2  f  2  ’  M’  T*  ^  d  a  ~  2 

aU  a 


(b- 15) 


and 
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(Tt  ->  N 1  ,  0 

tf  f  2  (a  ;  M,  -^)  d  a2  =  |  . 
O  rr 


(b- 16) 


2  2 

It  should  be  clear  that  both  and  o^T  are  functions  of  both  ^  and  M.  For 

large  values  of  M,  Holsinger  [  15]  provides  some  approximate  formulas  for 
the  confidence  interval 


2  _  2-1 

tfu  -  a  [1 


K 


-2 


JZ  (M-  1 ) 


] 


(b-  17) 


2  2 

crT  -  a  [  1  +• 


K 


-2 


V2(M-1) 


-] 


where  K  is  defined  by 


£_  - 
2 


uu 

-4zr  I 

JZn  K 


t2/2, 
e  7  dt. 


(b-  18) 


Thus  for  a  confidence  level  $/2  =  10$  and  for  a  ±  .2  dB  confidence  interval, 
the  required  number  of  samples,  M,  is  computed  from  equations  b-  17  and 
b-  18  to  be  about  1,  600.  As  normal  practice  we  have  chosen  a  sampling 
interval  of  50  msec  for  2  minutes  of  playback  data  so  that  approximately 
2400  samples  are  collected  for  the  variance  computation.  (It  should  be  noted 
that  50  msec  is  equal  to  five  time  constants  when  RC  =  ,01  sec.  The 
real  time  bandwidth  that  the  noise  is  subjected  to  is  Af1  =  100/F  =  3.4  Hz. 

The  wide  band  non-Gaussian  data  being  subjected  to  this  3.4  Hz  filter  comes 
out  reasonably  Gaussian.  ) 

An  example  of  the  output  of  the  noise  program  on  the  Varian  620i  is 
2  2 

shown  in  Fig.  B-2.  o'  is  given  as  1.265  (volt)  .  The  equivalent  input 

X 

noise  at  the  antenna  is  given  by  equation  b-9  which  is  rewritten  for  this 
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particular  experimental  configuration 


20  log  (H  ).T  =  [-128.91  +  20  log  ^  +  20  log  (j-)  +  20  log — 

0N  P1  fS  G 


+  20  log  — +  10  log  (CT  2)]  dBH/«/Hz, 

,2  X 


C 

s 

(b-  19) 


For  Utah  tape  no.  13  (at  1900  -  2000  GMT)  the  input  noise  intensity  at  45  Hz 
is  calculated  to  be  -134.2  dBH/  JTlz  ±  .  15  dB. 


io.oot 

I  1  -  6  -  1450?  | 

-  UTAH  TAPE  13 

^sompie  ^ 

RC  =  0.01  sec 


/?>  =  -6  dB 

TIME  INTERVAL  ( 1900  -  2000  GMT ) 
(  2  mm  R  B.  time) 

20  log  HN  ^  -  134.2  dBH/v/iHz 


-10  00  L 
0 


TIME 


2000 


Fig,  B-2.  Typical  noise. 
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APPENDIX  C 


Confidence  Interval  for  Signal  Estimates 

The  output  of  the  typical  receiver  structure  shown  in  Fig.  4  is  the 

A2 

random  variable  A  where 


A 


The  random  variables  ^(fi,  E )  or  (where  E  is  the  sample 

time)  appear  from  the  data  to  have  the  following  properties. 

(a)  At  each  sample  interval,  the  computed  variances  of  and  E 

c  s 

are  identical  within  experimental  tolerances; 

(b)  Since  the  equivalent  integration  bandwidths  are  of  the  order  of 

a  few  thousandths  of  a  Hertz  the  variables  E  and  E  are  very  close  to 

c  s 

being  Gaussian  distributed; 

(c)  The  means  are  A  cos  <p  and  A  sin  (p  respectively. 

Hence  the  probability  functions  are  therefore 


(E^  -  A  cos  <p 

2NT 


E 

c 


v'ZirNj 


fL  <Ls> 

S 


1 


(L  -  A  sin  <p)" 

Tn~t 


-v/2ffNT 


(c-2) 


The  next  assumption  made  is  that  and  Lg  are  uncorrelated.  Examina¬ 
tion  of  the  output  traces  of  the  individual  quadrature  channels  always  show  a 
lack  of  correlation.  It  is  suggested  that  the  extremely  narrow  filter 


bandwidth  is  the  cause  of  this  lack  of  correlation.  The  reasoning  is  as 

follows.  It  is  known  that  most  of  the  energy  is  contained  in  the  relatively 

infrequent  large  sferics.  These  sferics  pass  through  the  extremely  narrow 

filter  giving  rise  to  'transient1  sinusoid  with  a  certain  phase  relation  with 

respect  to  the  phase  of  the  receiver's  signal  reference  frequency.  Since  the 

phase  relationships  between  succeeding  large  sferics  and  the  receiver's 

reference  frequency  are  random;  one  could  expect  at  one  time  that  virtually 

all  the  noise  energy  shows  up  on  one  of  the  quadrature  outputs  whereas  the 

energy  from  a  succeeding  sferic  may  show  up  only  on  the  other  quadrature 

outputs.  A  randomization  of  the  phase  relationship  between  the  filtered 

sferic  and  the  receiver's  reference  frequency  clearly  leads  to  a  lack  of 

correlation  between  L  and  L  .  Hence  if  it  is  assumed  that  L  and  L  are 

c  s  c  s 

uncorrelated,  it  is  easy  to  show  (see  Van  Trees  [8])  that  the  cumulative 

^2 

distribution  function  for  A  is  related  to  the  Marcum  Q  function  by 

2  2 

a  oo  .  *.+  * 

P(A2sy2)=  1-J  UQ  (at)  e  2  d-t  (c-3) 

y  A/n7 

P(A2  ^  y2)  =  1  -  Q  ,  -^=)  (c-4) 

where 

l  = 

Let  us  select  a  confidence  level  ?  .  It  is  standard  practice  to  derive  the 
confidence  curves  from  the  relations 

A 

|  =  P(A2  sy2L)  (c-5) 
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and 

2  =  P<aZ  2  !&• 

Figure  C-l  shows  Marcum's  Q  function  with  the  curve  parameter  — — —  con- 
verted  to  dB.  Figure  C-2  shows  the  confidence  curves  (for  fixed  levels  of 
e  )  which  are  derived  from  equations  c-4  and  c-5  and  Fig.  C-  1.  Using  the 
confidence  curves  is  particularly  simple.  Assume  the  confidence  level  e  is 

/A 

specified  and  that  one  measures  a  value  A^/N,..  A  horizontal  line  is  drawn 

Z  Z 

on  Fig.  C-2  through  this  value  of  A  /N^.  The  two  values  of  A  /Nj  [defined 
by  the  intersection  of  the  horizontal  line  with  the  confidence  curves  (with 
parameter  $/Z]  form  the  limits  of  the  confidence  interval.  An  example  of 
how  to  construct  a  confidence  interval  is  also  shown  on  the  figure.  The 
confidence  interval  (A^/Nj,  A^/N^.)  itself  a  random  variable.  The  true 
value  of  A^/Nj.  is  included  within  this  interval  with  probability  1-^.  From 
these  confidence  curves  it  is  easy  to  see  that  the  confidence  interval  rapidly 
increases  in  size  with  decreasing  values  of  A  /N^. 

We  also  use  on  occasion  the  term  confidence  band.  To  determine  a 
confidence  band  one  first  computes  a  theoretical  signal  level  and  the  corres¬ 
ponding  signal  to  integrated  noise  ratio.  Then  a  vertical  line  is  drawn  on 
Fig.  C-2  through  the  value  of  A^/Nt  intersecting  the  appropriate  pair  of 
confidence  curves  at  two  points.  The  range  of  expected  sample  A^/N^.  values 
enclosed  by  two  horizontal  lines  through  the  two  intercept  points  will  be  called 
the  confidence  band.  In  brief  the  difference  between  the  confidence  band  and 
the  confidence  interval  is  the  following.  If  we  have  a  theoretically  calculated 
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Fig. 


Tl8-6-l3HI 


P(/>y)  *Q(a./9)  MARCUM'S  Q  FUNCTION 
A -PEAK  VALUE 


0  Z 


-  0.3 


I  2  3  4  5  6  7  8  9  10  H  12  43  M  (5  16 


C-l.  Cumulative  probability  distribution  (Marcum's  Q  function). 


1 14 


Pi/>Y)  = 


S/N  (<JB)  •  20  log 


Fig.  C-Z.  Confidence  curves  (for  Marcum's  Q  function). 
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signal  value,  then  the  confidence  band  gives  the  range  of  sample  biased 
signal  levels  that  one  might  experimentally  measure  {(1-^)  100$  of  the 
samples  would  be  in  the  band}  provided  that  the  theoretical  calculated  value 
was  in  fact  the  true  value.  If  on  the  other  hand  we  are  presented  with  a  set 
of  experimentally  measured  biased  signal  levels ,  then  for  each  of  these 
measurements  we  have  a  unique  confidence  interval.  The  statement  can  then 
be  made  that  with  probability  (1-^)  100$  each  interval  will  include  the  true 
value  of  signal  A  /Nj. 

Confidence  curves  (see  Fig.  C-2)  for  S/N  ratios  between  5  and  25  dB 

have  been  derived  from  the  Marcum  Q  function  tables.  However,  Marcum’s 

2 

Q  function  tables  are  not  tabulated  for  S /N  =  A  /Nj  >  25  dB.  Nevertheless, 

there  is  a  simple  asymptotic  expression  for  Q  ( — — —  ,  — * )  as  given  by  Van 

VNj 

Trees  [  1 6] 

(c-6) 

Q(-”~  ,  - j:  )  ^  f  1  +  er:f  (~ )  }  for  y  >  >  ^/Nj  and  y  >  >  y  -  A. 


Recalling  equation  c-4 


72Ni 


P(A2  £  y  )  =  1  -  Q(-  A 


y 


-v/Nj 


(c-7 ) 


it  follows  that 


P(A2  <;  y  )  =  j  fl  -erfA^)}  . 

V2NI 


(c-8) 


The  confidence  intervals  can  be  computed  from  this  distribution.  For 
example  assume  one  wants  to  determine  the  1-f  =  80$  confidence  curves. 
These  are  determined  from  setting  equation  c-8  equal  to  .  1  and  .  9  and 


determining  the  value  of 


A  -  y 


from  the  error  function  tables.  This  turns 


V2Nj 
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out  to  be 


<-£r- 


=  ±  1.286  for  1-f  =  80$. 


VNj  yNj 

In  table  C-  1  we  have  tabulated  typical  sets  of  confidence  curves  for  true 

S/N(=  A^/Nj)  ratios  between  23  and  35  dB.  It  is  notable  that  the  interval 

decreases  very  slowly  with  increasing  S/N  ratio. 

It  should  be  carefully  noted  in  this  section  that  A  is  a  peak  value  and 

not  the  rms  value.  Therefore,  since  the  computed  values  |H  (t.,f  )  L 

1  <p  i  s  'b,  1 

are  given  in  terms  of  rms  rather  than  peak  values  one  must  add  3  dB  to 
equation  19  before  the  confidence  curves  can  be  used  since  these  curves 
reflect  peak  rather  than  rms  values. 
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Table  C-l . 

Confidence  Curves  (1  -  e  =  80%)  for  S/N  greater  than  23  db. 


A/a  (db) 

A/a 

Lower  Limit 

Upper  Limit 

y/a  db 

(y/a  -  A/ a)  db 

y/a  db 

(y/a  -  A/a)  db 

23 

14. 125 

22. 17 

-.83 

23.75 

.75 

24 

15.849 

23.26 

-.74 

24.68 

.68 

25 

17.783 

24.35 

-.65 

25.61 

.61 

26 

1 9.  953 

25.42 

-.58 

26.54 

.54 

27 

22.387 

26.49 

-.51 

27.48 

.48 

28 

25.119 

27.54 

-.46 

28.42 

.43 

29 

28.184 

28.59 

-.41 

29.39 

.39 

30 

31.623 

29.64 

-.36 

30.35 

.35 

31 

35.482 

32 

39.812 

31.71 

-.29 

32.28 

.28 

33 

44.669 

34 

50.120 

35 

56.236 

34.80 

-.20 

35.20 

_ 

A  is  the  true  value  of  the  signal  level 
a2  1/2 

(A  )  '  is  the  sample  value  of  the  signal 
a  =  J N |  is  the  square  root  of  the  integrated  noise. 
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APPENDIX  D 


Derivation  of  Composite  Confidence  Limits  for  A^(E)  and 

In  Section  III  we  defined  the  unbiased  composite  field  strength  at  the 


frequency  f.  as 


£  (A2(fi))u.i 


<a2<v>  =A 

S  (t 


N 


LL 


(d-1) 


•  1  nt  . 

1=1  i,  i 


ys 

2, 


where  N  .  is  the  ith  integrated  noise  sample  estimate  and  (A  (f.))  .  is  the 

■!»  1  J  u»  1 

corresponding  ith  unbiased  estimate  for  A  2 ( f. ) .  It  is  noted  that  the  weighting 
gives  greater  weight  to  those  samples  which  exhibit  the  higher  S/N  ratio. 

The  Nt  .  can  be  determined  to  a  very  high  degree  of  accuracy  and  for 
purposes  of  computing  a  composite  confidence  limit  these  values  will  be  con¬ 


sidered  as  known  constants.  The  mean  and  variance  of  A  (f. )  then  are 

J 

M  A 

£  E[A2(f.)l 

-  •  ,  Nt  .  L  j  JU,1 

i=l  I,  i 


E[ A  (f. )]  = 


M  1  2 

?  TT" A  (fj: 

i=  1  I.i  J 


(d-2) 


1 


M 

T  i  nt  . 

1=1  i,  i 


i 


M 

1  NT  ' 
1=1  I,  1 


■ A  <V 


and 


M  ^ 

,  ?  .(NTT?  var  iA 

2,,  vx  _  i=l  I,  i _ J 


var(A  (f. ))  = 


M 

(£ 


1 


•  i  Nr  • 
i=l  I,  l 


(d-3) 


2  2 

It  is  seen  that  A  (f.)  is  an  unbiased  estimator  for  A  (f. ).  For  the  receiver 

J  J 

structure  shown  in  Fig.  4  it  can  be  shown  that 
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var  [A2(f.)l  .  =  var  [A2(f.)l,  .  =  4(N2  .  +  NT  .  A2(f.))  . 

1  v  J  Ju,i  1  J  Jb,i  I,i  I.i  J 


For  A^(fj)/Nj  j  >  10  we  can  make  the  approximation 
xx 

var  [A2(f.)l  •  =“  4A2(f.)NT  . 

1  J  J  u,  i  J  I.i 


(d-4) 


and  hence 


M  4A  (f.) 
L  - L 


var  (A*(f.)) 


•  -I  NT  . 

i=  1  I,  l 


j-  M 

N~ 

1=1  1,1 


4A2  (f 


M 


1 -  £ 


1 


(L  1  nt  •  ' 

1=1  1,1 


=  k.  A“(f.: 

J  J 


(d-5) 


where  K 


J  M  i 


T  >T  ' 

1=1  I,  1 


Note  that  the  noise  Nj.  ^  is  the  ith  sample  of  integrated  noise  centered  at  the 
frequency  F.  The  random  variable  A^(F)  can  be  assumed  to  be  Gaussian 
distributed  in  the  vicinity  of  the  mean  A^(F).  In  equations  d-6  through 
d-12  we  assume  that  the  dependence  on  F  is  implicit  so  as  to  simplify  the 
notation.  Hence  the  probability  density  function  is 


f  (A2  ;A2 )  = 


ka 


(A 2  -  A2)2 

2  (AK)2 


(d-6) 


The  confidence  curves  are  then  derived  from  the  relationship 


(A2  )jj  =  y+ A2 


L 


(A2)t  =  -y+A2 


f(A2  ;A2)dA2  =  1  -  e  . 


(d-7) 


A2  -  A2 

Let  us  normalize  by  setting  y  =  — -T-gr -  .  Equation  d-7  reduces  to 
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Y 


/AK 


I 


-  y  /ak 


[  e  Y  '  2]  dy  =  1  -  c 

J2n 


The  upper  and  lower  bounds  are  both  determined  from  the  following 


Q<^>  = 


±  00 


J 


Y 


/AK 


{  ZZT 

V2ir 


e"y2/2}  = 


r/2 


(d-8) 


7/2  =  . 30 

f/2  =  •  10 

o 

ii 

f/z  =  .  025 

Y  /AK=  ft  ^ 

B  =  .524 
^  € 

ft  =  1.28 

H  c 

=  1.64 

B  =1.96 
c 

The  upper  and  lower  confidence  curves  are  thus  derived  from 


(A2  )y  =  A2  +  ft  (AK) 

(A2)  =  A2  -  ft  (AK) 

-b'  f 


(d-9) 


Figure  D- 1  shows  a  typical  composite  confidence  curve  for  the  averaged 
daytime  and  nighttime  unbiased  signal  estimates  at  f  =  45  Hz.  The  corres¬ 
ponding  confidence  intervals  (at  ^/Z  =  .  1)  for  the  composite  values  (A^(45))^ 
_ 

and  (A  (45))^  are  also  shown. 

The  last  set  of  confidence  curves  one  needs  to  compute  are  those  for  the 
2 

random  variable  A  where  we  have  defined  (f  is  the  mid- point  or  center 

m  r 

frequency)  _ 


A2(f  ) 
'  m 


1 


m  '2J+  1 


2  J  +  1 
)  L 
j=l 


A2  (f  • ) 

J 

f2 

J 
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Fig.  D-l.  Composite  confidence  curve. 
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Now 


=f=  2  ,  2J+1  E[A*(f  )] 

ElA‘ (fm>l  =  4  <2J+I>  E,  f2  1 

J=1  'j 


Define 


E(a2(V1  =  A2t6A' 


(d-  10) 

2 

A  is  independent  of 


2  2 

where  6  A.  is  a  small  perturbation  compared  to  A  . 

the  frequency  f.  and  related  to  the  direct  component  |H^(fm)  |^.  (Recall 

2  2 
that  A  jf. )  is  related  to  the  total  field  component  |H  (f.)  I  and  varies  with 

J  <p  J  1 

one's  position  with  respect  to  the  standing  wave  pattern.  )  Equation  d-9  can 

be  written  as 


2J+1 


2J+1  6  A 


E[A2(f  )1=  f2  (t-FTT)  {£  4"  +  S  — 5— } 

1  '  m  J  m  2J+1  l. ,  ,2  .  ,  .2  J 


j  =  1  f  7  j  =  1 

J  J  J 


f2 

J 


(d-11) 


For  the  receiver  locations  deployed,  the  data  show  that  in  general, 

2  2 

10  log  (A  /6Aj)  >  >  10  dB.  Moreover  since  we  have  designed  the  experiment 

to  obtain  (2J+1)  equispaced  samples  over  one  standing  wave  length  one  can 

2 

be  reasonably  sure  that  the  6  A^  fs  are  fairly  evenly  distributed  about  zero 

and  that  the  first  summation  in  equation  d-11  swamps  the  second  and  hence 

z  2  T4- 1 

2  2  1  x  1  ?  A  2 

E1A  <fm»  ~  fm  A  =  CA  ' 

J=‘  fj 


Since  C 


1  (i.e.  ,  C 


1.01  for  the  45  Hz  band)  we  can  set  C  =  1  and  hence 

(d- 12) 


2  2 

E[A  (f  )]  =*  A  or  equivalently 


E[  [H2 (f  )]  -  |H2(f  )  I  . 

1  I  <pv  m'J  I  <pv  m'  If) 
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The  last  equation  states  that  (f  )  |  is  an  unbiased  estimator  for  the  direct 


O 


field  component  |H  (fm)  |^. 


We  can  also  compute  the  variance  of  A  (f  )  since 


and 


-=77  4  ,  ?  2  J+  1  2J+1  E[A2(f.)A2(f .)] 

e'a  »  =  4 'ik>  ,s. ,  =  ,  —&r — 

1  J  i  j 


var  (a2)  =  fm  E 


2J+ 1  2J+1 


i=l  j=l 


E[A2(£.)A2(£.)]-E[A  (£ .)]  E  [  A2  (f  -  >] 

«  f — — 5 — rh - 1 - ■- 

f  f 
1  J 


}  * 


A  (f.)  and  A  (T)  (where  i  /  j)  are  constructed  with  signal  estimates  taken 
from  two  completely  different  set  of  days.  Hence  there  is  a  lack  of 


correlation  in  the  noise  and 


E[A2(f.)A2(f.)]  =  E[A2(f.)]  E[A2(f.)] 


for  i  /  j. 


This  results  in 


— 7  4  1?  2J+1  var  [A  (f.)} 

var<A  >  =  4  <2JFT>  =  —4 - ^ 

t 


(d-  13) 


However  from  equation  d-5  we  have 


var  {A^(fj)}  =  Kj^A^(T)  where  j  refers  to  the  jth 


frequency  and  hence  in  the  same  way  as  above  one  can  make  the  approximation 

/  I  4.  I  K  U  “  IT  ■  ✓  I  _l_  I 

2,  A 


var(A  )  =  fm  2 

J=  1 

=  k2a2 


,  ,  2J+1  K2A2(f.)  .  ,  ,  2J+1  Kr 

1  ,2  ^  J - L2_^f4(  1  }2(  — J- )  A 2 

m  '2  J+ 1  ,  r4  ’ 

1=  1  f. 

J  J 


124 


where 


2  A  A 
K  =  f 

m 


1 


2J+1  K. 

2  1  E  "4"}  • 

(u+if  j=i  y 


(d- 14) 


Again,  assuming  that  A  is  Gaussian  distributed  also  one  can  use  as  the 

2  2 

distribution  function  equation  d-6  with  (A,K  )  being  replaced  by  (A,K  ). 
Then  the  composite  confidence  limit  for  A^  can  be  constructed  in  the  same 


wa  y  a  s  tho  s  e  f o  r  A  ( f  j  ) . 
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APPENDIX  E 


Receiver  Calibration 

A  necessary  condition  for  accurate  estimation  of  a  is  that  each 
receiver  system  be  well  calibrated.  We  are  looking  for  estimates  of  a 
accurate  to  within  ±0.1  dB/Mm  so  it  is  necessary  that  each  station’s 
receiver  be  calibrated  to  less  than  ±  0.  3  dB.  That  is  in  the  absence  of 

noise  one  must  be  able  to  estimate  the  signal  H  to  within  that  accuracy. 

0 

It  is  apparent  that  the  errors  in  meter  readings,  in  taking  values  from  the 
chart  recordings,  in  calibrating  the  chart  recorders  and  PAR’S,  etc.  ,  when 
summed  together  cannot  exceed  ±0.3  dB.  This  requires  an  extremely 
exacting  calibration  routine  which  is  described  in  the  remainder  of  this 
section.  As  noted  in  the  text  each  complete  receiver  is  composed  of  two 
parts;  the  field  site  recording  equipment,  and  the  laboratory  based  play¬ 
back  facility.  The  laboratory  based  playback  facility  is  of  course  common 
to  all  the  receivers.  The  remainder  of  this  appendix  first  outlines  the  field 
site  equipment  calibration  and  then  discusses  the  calibration  of  the  laboratory 
playback  facility. 

Calibration  of  the  Field  Site  Receiver  Segment 

The  primary  function  of  the  calibration  is  to  determine  the  system 
gain  at  the  received  signal  frequency  and  to  determine  any  changes  in  this 
gain  over  the  recording  duration.  For  this  purpose  a  calibration  signal  is 
continuously  injected  in  series  with  the  antenna  (refer  to  Figs.  4  and  5  in 
Section  III).  The  calibration  frequency  is  selected  to  be  1  Hz  removed  from 
the  received  signal  frequency.  Its  amplitude  was  chosen  to  be  sufficiently 
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large  such  that  its  integrated  signal  to  noise  ratio  exceeds  30  dB.  The 
following  notation  is  used;  G(f)  is  the  open  circuit  voltage  of  the  loop  antenna, 
B  (f)  is  the  calibration  voltage  impressed  across  the  0.  5p  resistor,  A  is  the 
signal  amplitude  input  into  the  quadrature  correlators  in  the  signal  channel, 

B  is  the  calibration  amplitude  input  into  the  quadrature  correlators  of  the 
calibration  channel  and  VQ(f)  is  the  voltage  output  of  the  narrowband  channel 
before  the  tape  recorders.  Let  G  represent  the  voltage  gain  from  the 
transformer  input  to  the  output  of  the  narrowband  filter.  Thus  with  the 
antenna  short  circuited  and  by  successively  injecting  into  the  calibration 
pad, sinusoids  of  frequency  fg  and  f^,  one  can  then  measure  the  gains 


A  Vo(fG: 


e  (fc) 
VJfc 


r  A  _ 

s  -  B  (f 


o'  S 


ant  S.  C. 


ant  S.  C. 


-10, 


where  B  (f)  =  Pj.VCM(f)  =  (4.4959  X  10  )(VCM(f)[mv.rms] 


(e-1) 

(e-2) 

(e  -  3 ) 


is  the  voltage  attenuation  ratio  of  the  precision  attenuator  pad.  Both 
^CM  anc*  Vo(f)  are  measured  with  a  precision  rms  voltmeter.  The  gain  of 
the  combined  record  circuit  of  the  P.  I.  tape  recorder  and  the  playback 
circuit  of  the  Ampex  tape  recorder  was  measured  and  found  to  be  independent 
of  frequency  (i.  e.  ,  G^(f^)  =  G^(fg)  =  G^).  The  system  gains  with  the 
antenna  included  in  the  circuit  are  denoted  by  G^  and  Gg.  The  overall 
system  gain  from  the  antenna  (or  attenuator)  to  the  input  to  the  correlator 
receiver  is  given  by 


1Z7 


(e-4) 


A(fq) 

_ 2  —  s"'  t  p 

G  ( f c )  ~  US-UT 


B(fr) 

- —  =  G 1  G 

e(f~)  C  •  • 


However  since  the  antenna  impedance  is  small  compared  to  the  transformer's 
input  impedance 


GC  = 

Gs  ’  Gs 


(e-5) 


Equations  e-4  and  e-5  are  then  combined  to  form  the  ratio 


G(fs) 


A (f  )  G 

-  ( - —  )  / - 

«(IC)  ,B(Ic)"gs)- 


(e-6) 


Recall  from  equation  a-2  that  the  input  field  threading  the  antenna  loop 
is  related  to  the  open  circuit  voltage  G(fQ)  by 


a  (fs) 

=  (.  007  17 8 )fc 


(e-7) 


With  noise  present,  the  outputs  of  the  correlation  receivers  are  the  estimates 

A  /s 

z  z  . 

A  and  B  .  Then  equations  e-3,  e-6  and  e-7  can  be  combined  to  yield  an 


equivalent  input  H  field  estimate 


rn 


J H7 


<P 


P"Wc) 

S  B 


1 

(.  007  178 )fc 


The  following  alternative  form  is  more  useful 

A  G 

10  log  H2  =  -144.  06  +  20  log  -  20  1og  fs  +  201og  VCM(fc) 

1  S 

+  10  log4j  (e  - 8 ) 

B 
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where  V-w  is  in  mv.  (rms). 

CM 

In  summary,  the  field  site  calibration  consists  of  measuring  and 

GQ  via  equations  e-1,  2  and  3  both  at  the  beginning  and  at  the  end  recording 

o 

period.  V  .,(f  )  is  also  measured  at  these  times.  Average  values  (  ) 

JVl  V_> 

of  these  three  quantities  are  taken  if  there  are  small  changes  in  value  over 

2 

the  recording  period.  It  is  noted  that  B  is  continuously  plotted  on  the  chart 
recorder  of  the  laboratory  playback  facility.  This  provides  a  continuous 

check  of  the  system  gain  during  the  recording  period.  The  four  quantities 
~  ~  ~  2 

(G^,  Gg,  V£^.(f(g)  and  B  )  are  the  calibration  constants  to  be  used  in 
equation  e-8  for  the  field  strength. 

Calibration  of  the  Laboratory  Segment  of  the  Receiver 

The  main  elements  to  be  calibrated  in  the  laboratory  playback  segment 
of  the  receiver  are  the  correlation  averagers  (i.e.  ,  the  Princeton  Applied 
Research,  Model  120  lock-in  amplifiers ) ,  squarers,  summers  and  chart 
recorder.  The  main  difficulty  is  calibration  of  the  PAR's  and  in  particular 
maintaining  each  pair  of  PAR's  in  exact  quadrature.  One  finds  that  if  one 
relies  on  the  stepped  phase  switch  (0°,  90°,  180°  and  270°)  on  the  PAR's  a 
possible  error  of  the  order  of  1  dB  in  signal  estimate  can  result.  (This 
error  is  a  consequence  of  the  fact  that  the  phase  steps  are  not  as  close  to 
90°  as  advertised.  )  Figure  E-1  shows  an  operational  equivalent  of  a  single 
PAR  and  Fig.  E-2  shows  the  calibration  set  up  for  a  pair  of  PAR's.  For 
frequency  and  input  scale  selected  the  following  procedure  is  followed 

a)  With  the  calibration  oscillator  set  at  the  desired  frequency  and 

amplitude  (such  that  the  reading  on  the  precision  rms  voltmeter  corresponds 
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Fig.  E-l.  PAR  representation. 
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Calibration  network 


to  full  scale  input  voltage)  adjust  the  phase  shifter  so  that  the  phase  meter 
reads  0°.  On  PAR  no.  1  adjust  first  the  frequency  trim  controls  of  the  tuned 
amplifiers  so  that  the  output  reading  on  the  digital  voltmeter  (DVM)  is  a 
maximum.  Then  adjust  the  gain  adjust  control  so  that  the  output  reads  a  full 
10.0  volts.  Ignore  PAR  no.  2  for  the  present. 

b)  Adjust  the  phase  shifter  so  that  the  phase  meter  reads  90°.  On 
PAR  no.  2  perform  the  same  adjustments  as  performed  on  PAR  no.  1  in  the 
preceeding  step.  PAR  no.  1  should  have  a  reading  close  to  0  volts.  Adjust 
the  variable  phase  on  PAR  no.  1  a  few  degrees  so  that  the  output  does  in 
fact  read  exactly  zero. 

c)  Adjust  the  phase  shifter  so  that  the  phase  meter  again  reads  0°. 
The  output  of  PAR  no.  1  should  read  10V.  If  not  a  small  gain  control  adjust¬ 
ment  should  make  it  read  10V,  PAR  no.  2  should  read  close  to  0  volts. 

Adjust  the  variable  phase  control  on  PAR  no.  2  so  that  it  does  read  0  volts. 
The  two  PAR'S  now  should  be  individually  calibrated  and  in  quadrature. 

d)  As  a  final  check  with  the  DVM  connected  to  the  output  of  the 
summer  sweep  the  phase  shifter  through  360°  noting  the  output.  If  properly 
calibrated  the  output  should  not  vary  by  more  than  a  0.  1  dB. 

The  squarers  and  summer  can  be  easily  adjusted  so  that  each  of  the 
segments  agree  with  the  predicted  behavior  to  within  0.  1  dB.  The  chart 
recorder  is  for  each  data  run  calibrated  indirectly.  That  is  after  each  run, 
with  zero  inputs  to  the  PAR's,  on  ore  PAR  a  range  of  bias  voltages  (Measured 
with  the  DVM)  are  selected.  The  corresponding  values  then  are  marked 
on  the  chart  recording  next  to  the  point  to  which  the  recorder  pen  was 
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deflected.  Hence  on  the  chart  we  will  have  several  values  for  the  scale  of 
the  amplitude  A  (not  the  scale). 

Another  practical  point  concerns  the  behavior  of  the  PAR'S,  Tape 
speed  variations,  if  not  compensated  for,  can  yield  an  error  in  signal 
estimate  of  1  dB  or  greater.  This  is  due  to  the  fact  that  the  frequency 
responses  of  the  tuned  amplifiers  in  the  signal  PAR  pair  and  in  the  calibra¬ 
tion  PAR  pair  are  not  identical.  Speed  variations  of  the  order  of  5^  can  give 
errors  of  the  order  of  1  to  2  dB  depending  on  the  particular  band  pass  filters 
used  in  the  tuned  amplifiers.  This  problem  has  been  circumvented  by 
addition  of  a  speed  control  loop  which  compensates  within  1$  for  speed 
variations  in  both  the  recording  and  playback  tape  recorders. 
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APPENDIX  F 


Estimates  of  the  Effects  of  the  Round- the -World  Wave 

It  is  apparent  that  the  earth-ionosphere  cavity  can  have  a  standing  wave 
pattern  with  significant  nulls  and  peaks,  wherever  the  direct  and  round-the- 
world  signal  components  are  comparable  in  magnitude.  This  is  normally 
the  case  over  a  wide  range  of  distances  in  the  vicinity  of  the  antipode  for 
propagation  in  the  lower  ELF  domain.  Assuming  a  world-wide  isotropic 
and  homogeneous  ionosphere  of  constant  height,  the  total  field  is  simply  the 
sum  of  a  direct  component  and  an  around-the- world  component  given  by 


(H  )_  =  (H  )  +  (H  ) 


<P 


(pi 9  cp  ft 


where 


jk  d 

,,  J  o'v  , 

(H  )  =  F  e~a  d  e  ph 

<n  & 


(f-D 


and 

-a’(2 nre-d)  jkQ(^— )  (Z^^-d) 

(H  )  =  -F  e  e  ph  .  (f-2) 

O  TC 

F  is  given  by  equation  3  (Section  II)  and  d  <  Trr^.  Since  the  ionosphere  is 
not  isotropic  and  homogeneous,  then  a  better  approximation  is  to  replace  the 
exponential  dependence  by 
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'd  ->  f  a’£( r)dr  ; 


ph 
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dr 


JZ 7Tr  -d) 

'(2nr  -d)  -  J 


Of 1  ( r )d r  ;  etc. 
W 


resulting  in  equations  6  and  7  in  Section  II.  It  must  be  mentioned  that  at 
present  there  is  no  complete  mathematical  justification  for  the  simple 
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replacement  of  exponential  terms  (although  it  is  consistent  with  one's 

physical  intuition).  However,  it  tends  to  accentuate  the  nulls  and  peaks  of 

the  standing  wave  pattern  and  hence  provides  good  bounds  for  the  true  value 

(H 

of  total  signal  strength  (H  )rr .  Figure  F-l  shows  curves  of  I — ^ — I  and 

(HI  .  9 

I — ff  rcu  for  various  constant  values  of  a  and  a  .  In  absence  of  phase 
1  F  1  -  R  ^ 

information  one  can  compute  the  maximum  and  minimum  values  of  from 
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Figure  F-2  is  a  graph  of 
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versus 
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We  use  these  curves  in  the  following  ways: 


(1)  If  the  attenuation  coefficients  a  and  a  are  known  constants 

&  R 

and  do  not  depend  on  location  one  can  use  curves  similar  to  those  in 

Fig.  F-l  to  determine  |(H  )  A  l/F  and  l(H  )  l/F  as  well  as  the  ratio 
| (H  )  |  '  '  “>»' 

I'/TT  \ ^  1 •  If  these  are  known  then  the  maximum  and  minimum  values  of  (H  )m 

I'V*1  ®T 

can  be  found  from  Fig.  F-2. 

(2)  If  the  attenuation  coefficients  (v  and  a  vary  with  direction 

JS  ft 

and  position  then  one  must  use  the  full  equations  6  and  7  of  Section  II.  Then 

with  the  use  of  Fig.  F-2  one  can  compute  the  maximum  and  minimum  values 

of  (H  ) 

<P  T 
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Fig.  F-l. 


Normalized  direct  and  round 


the  world  field  strength 
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Fig.  F-2.  Maximum  and  minimum  standing  wave  interference  effects. 
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APPENDIX  G 


Location  of  Standing  Wave  Minima  and  Maxima  and  Magnitude  of  Standing 

Wavelength 

Since  we  propose  to  change  the  position  of  the  standing  wave  pattern 
relative  to  the  receiver  stations  by  changing  frequency  in  small  intervals  it 
is  necessary  that  we  obtain  reasonable  assurance  that  the  pattern  shifts  by 
about  a  full  standing  wave  length  over  the  frequency  range.  If  this  is  so, 
the  average  of  the  total  field  strengths  (weighted  appropriately)  will  be  a  good 
approximation  to  the  direct  field  strength  at  the  center  frequency  (in  this 
case  at  45  Hz ) . 

If  one  chooses  as  a  model  for  the  earth-ionosphere  cavity  the 
homogeneous  cavity  considered  in  the  first  part  of  Section  II  it  can  be  shown 
that  the  minimums  of  the  |H^  |  ^  standing  wave  pattern  due  to  a  surface 
dipole  radiator  are  displaced  from  the  antipode  by  a  distance  p1  where 

n  =  0,  1 ,  Z , .  .  . 

The  standing  wave  length  is  then  found  from  the  above  expression  to  be 


(4n  +  3)c 

p  = 


*sw  =  P* 


n+l  "  pn  2(c/vph)  2f(C/vph) 


Now  Hawaii  is  the  station  furthest  away  from  the  transmitter  (d  =  6.  53  Mm 
and  13.  47  M  away  from  the  antipode)  and  hence  the  station  most  effected  by 
the  interference.  As  an  example  let  us  look  at  the  position  of  fourth 
minimum  from  the  antipode  (for  c/v^  =  1.Z5  and  f  =  45  Hz) 

g 

p'( 45,  n=4)  =  8(45)1  26 — ^  =  56  Mm  from  the  antipode  and 
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X  =  2.  64  Mm. 
s  w 

Now  assume  we  change  the  frequency  to  46  Hz  and  again  assume 
0'(46,  n=  4)  =  12.27  Mm  and 


)  = 


1.26 


X  =  2.  58  Mm. 
s  w 

Thus  although  the  standing  wavelength  changed  only  by  .  06  Mm  the  position 
of  the  standing  wave  pattern  changed  by  .30  Mm  which  is  greater  than  10^  of 
the  standing  wave  length.  At  distances  nearer  to  the  transmitter,  the  shift 
in  the  null  will  have  even  a  larger  magnitude.  It  should  now  be  obvious  that 
with  a  relatively  small  change  in  frequency  and  hence  in  standing  wave  length 
we  nevertheless  get  a  significant  change  in  the  position  of  the  standing  wave 
null.  Due  to  the  structure  of  the  transmitter’s  modulator,  it  is  constrained 
to  change  frequency  in  steps  of  1  Hz.  Hence  the  next  step  is  to  determine 
the  number  of  frequencies  that  we  will  require  to  make  a  reasonable 
approximation  for  |H  at  45  Hz. 

Figure  G-l  is  a  general  set  of  curves  which  show  how  the  maxima  move 
both  as  a  function  of  frequency  and  c/v  The  minima  (which  are  not 
indicated)  are  half  way  between  the  maxima  curves.  The  vertical  axis  is  the 
distance  from  the  transmitter  antipode  normalized  by  the  free  space  wave¬ 
length  X  =  c/f.  It  should  be  evident  by  looking  at  this  figure  that  for  any 
fixed  value  of  c/v  if  we  change  frequency  over  a  10  Hz  bandwidth  between 
40  and  50  Hz  the  position  of  the  standing  wave  pattern  will  change  by  more 
than  a  standing  wave  length  with  respect  to  either  of  the  site  locations.  (In 
fact  we  could  even  select  only  those  integer  frequencies  41  to  49  Hz  and  be 
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NORMALIZED  DISTANCE  FROM  RECEIVER  TO  THE  TRANSMITTER’S  ANTIPODE  (p'/Z) 


Fig.  G-l.  Location  of  standing  wave  maxima. 
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reasonably  sure  that  we  have  a  shift  of  more  than  a  full  standing  wave  length.  ) 
The  argument  can  even  tolerate  small  changes  of  (c/v^)  with  frequency. 

Accordingly,  we  have  selected  a  five  frequency  experiment  with 
frequencies  at  41,  43,  45,  47  and  49  Hz.  The  frequency  band  then  is  8  Hz 
(~  18$  center  frequency  of  45  Hz).  The  expected  decrease  of  c/v  ^  over 
this  band  is  of  the  order  of  only  2$.  Hence  there  is  little  risk  in  equation  4 
that  the  product  ^c/v^)  is  a  constant  over  the  frequency  range  and  as  a 
result  the  step  frequency  changes  should  result  in  significant  shifts  of  the 
standing  wave  pattern. 

The  preceding  comments  can  also  be  applied  to  the  7  0-80  Hz  range  in 
an  attempt  to  determine  say  |H^(7  5)|^.  The  number  of  transmissions  in  the 
70-80  Hz  range  was  only  1/3  of  those  in  the  40-50  Hz  range.  As  a  con¬ 
sequence  no  attempt  was  made  to  perform  as  comprehensive  an  analysis. 
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APPENDIX  H 


Transmission  Schedule 

Due  to  matters  of  economy  the  transmitter  was  operated  only  for  two 
eight  hour  periods  per  day.  One  period  was  chosen  to  maximize  the  number 
of  daylight  hours  for  the  entire  Noval  Scotia  to  Hawaii  path  and  the  other 
period  the  number  of  nighttime  hours  for  the  same  path.  We  thus  have  pur¬ 
posefully  avoided  whenever  possible  the  sunrise  and  sunset  periods  (where 
we  know  the  signal  level  changes  rather  rapidly).  The  Nautical  Almanac 
(1971)  can  be  used  to  give  local  sunrise  and  sunset  times  for  the  period  of 
experimentation.  (An  equally  useful  set  of  tables  which  also  include  the 
high  latitudes  is  reference  [  17  ]  .  )  The  time  which  we  consider  ionospheric 
sunrise  or  sunset  corresponds  to  when  the  sun's  zenith  angle  X  =  99°  (this 
would  correspond  to  the  time  when  an  observer  who  is  stationed  at  an  alti¬ 
tude  of  85  km  would  begin  to  observe  the  limb  of  the  sun  over  the  bulge  in 
the  earth).  The  estimated  ionospheric  sunset  and  sunrise  times  are  roughly 
interpolated  from  the  local  astronomical  twilight  and  ground  sunset  (sunrise) 
times.  Longitudinal  corrections  are  made  and  then  all  the  times  are  con¬ 
verted  to  GMT.  Table  H-l  shows  the  times  for  the  reference  date  March  26, 
1971. 

Table  H-2  shows  graphically  the  day-night  segments  for  all  the  sites 
and  hence  gives  an  indication  of  how  the  transmission  times  were  chosen. 

An  effective  integration  time  T  =  1.  42  hrs.  was  chosen  for  the  bulk  of  this 
experiment.  The  position  of  the  corresponding  sample  times  with  respect 
to  the  relevant  all  day  or  all  night  path  is  shown  as  well. 
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The  transmission  intervals  were  thus  scheduled  for 


0300  -  1100  GMT  (Nighttime  transmissions) 
1615  -  0015  (Daytime  transmissions) 


Toble  H-1.  Sunrise  ond  Sunset  Times 


Ref.  Dote  March  26,  1971 


Sunrise  Sunset 


Location 

Lotitude  °N 
Longitude  °W 

Begin 

Astronomicol 
Twi  1 ight 

Est. 

Ionospheric 

Sunrise 

Sunrise 

X  =  91° 

Sunset 

X  -  91° 

Est. 

Ionospheric 

Sunset 

End 

Astronomicol 
Twi  1  ight 

Howoii 

19°  15'  N 

0447 

0523 

0600 

1813 

1850 

1927 

155°  15'  W 

0508  f 

0544  f 

0621 

1834 

1911 

1948 

1544  GHT  %/ 

0511 

Wendover 

40.66° 

0425 

0510 

0555 

1817 

1857 

1947 

Utoh 

114.1° 

0401 

0446 

0531 

1753 

1833 

1925 

1246 

02  33 

p-Sife 

46°  04' 

0415 

0504 

0553 

1819 

1909 

1959 

Wisconsin 

90°  53' 

0419 

0508 

0557 

1823 

1913 

2003 

1 108 

0113 

Nova  Scotio 

44°  26' 

0415 

0504 

0553 

1819 

1909 

1959 

65°  13' 

0429 

0518 

0607 

1833 

1923 

2013 

0918 

2323 

X  =  Zenith  Angle  of  Sun 

Astronomical  Twilight  ot  X  =  108° 

Est.  Ionospheric  Sunset  or  Sunrise  X  -99° 
(corresponds  to  h  =  85  km) 


*  Sfd  Locol  Time  for  Specified  Lotitude  ond 
Longitude  of  15°  (n)  n  -  0, 1 , . . .  23 

t  Sfd  Locol  Time  with  Longitude  Correction 

v/GMT  with  Longitude  Correction 


*Volues  ore  from  Americon  Epherimis 
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TIME  (GMT) 


TABLE  H-2 


0000 


0200 


0400 


0600 


0800 


1000 


1200 


MOO 


1600 


1800 


2000 


2200 


2400 


Transmission  Periods 


Ionospheric  Parameters 
an  26  Mar  1971 


Sample  Times 


Nova  Scatia  -  Utah  Path  Utah -Hawaii  Path 


-£-0015 


-0300 


~  £ 


“0113 


t0918 


- 1100 


J0233 


t1108 


(/) 

00. 


2323 


0511 


E  I 
£  I 


I! 


8  I 


-0450 


—  0615 


—  0748 


-0918 


-I-  -1246 


1544 


-  1805 


-  1930 


-2115 


-2300 


I 


IT 

I! 

z  I 


-0615 


|  -0748 


sj 

t/>  I 

*  I  _ 


$1 

£4. 


0918 


1100 


-  1544 

-  1805 

-  1930 

-2115 

-2330 
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APPENDIX  I 


Error  Budget 

The  following  is  a  discussion  of  the  known  error  sources  which  con¬ 
tribute  to  error  in  a  single  measurement  at  a  particular  site.  It  should 
be  noted  that  signal  estimates  are  usually  based  on  a  number  of  individual 
measurements,  the  effect  of  which  is  to  average  some  of  the  errors  dis¬ 
cussed  herein.  This  effect  will  be  discussed  later. 

Turns  Area  Product  -  The  number  of  turns  on  the  antenna  is  known 
exactly  (471).  It  is  estimated  that  the  mean  diameter  can  be  measured  to 
within  ±.  25  inches.  Estimated  error  due  to  inexact  turns  area  product  is 

0.  8g. 

Calibration  Pad  -  An  attenuation  pad  is  employed  to  reduce  the  in¬ 
jected  calibration  signal  to  a  known  level  and  is  injected  in  series  with  the 
antenna.  The  attenuation  is  constructed  with  a  combination  of  .  01  and  .  05^ 
resistors.  The  overall  estimated  error  is  0.  1^. 

Calibration  Level  -  Field  Measurement  -  The  calibration  level  applied 
to  the  calibration  attenuator  is  measured  in  the  field  and  is  used  to  determine 
the  system  gain.  As  discussed  elsewhere  the  meters  are  calibrated  against 
a  secondary  standard  at  the  voltage  where  they  will  be  employed.  In  view 
of  this,  the  maximum  error  will  accrue  from  meter  reading  and  is  estimated 
as  1.0^. 

Frequency  Gain  Difference  -  The  calibration  signal  is  injected  1  Hz 
removed  from  the  received  signal  frequency  and  used  to  calculate  system 
gain.  Since  the  gain  at  the  two  frequencies  is  not  necessarily  the  same  a 
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correction  is  required  to  properly  reference  the  two  gains.  This  measure¬ 
ment  is  made  with  a  shorted  antenna  at  the  field  site  before  and  after  a 
measurement  period.  Calculations  are  based  on  the  average  of  the  two 
measurements.  Error  is  estimated  as  2^. 

Calibration  Level  -  Laboratory  Measurement  -  The  calibration  signal 
must  be  estimated  from  the  tapes  during  playback  at  the  laboratory.  Two 
factors  enter  into  the  error  calculations:  (1)  Calibration  of  the  instrument 
which  is  used  to  estimate  the  level  and  (2)  accuracy  of  the  squaring  circuits 
used  to  form  sum  of  squares.  A  signal  substitution  method  is  employed  to 
estimate  the  calibration  level.  Under  these  conditions  a  calibrated  voltmeter 
is  used  to  measure  the  substitution  signal  and  the  error  is  due  to  meter 
reading  (1^).  The  squaring  circuits  have  an  accuracy  of  1^  and  two  are 
employed.  Assuming,  worst  case,  that  all  errors  are  cumulative  the 
error  due  to  this  effect  is  3^. 

Signal  Level  -  Laboratory  Measurement  -  The  same  methods  are  em¬ 
ployed  in  signal  estimation  as  for  calibration  level  and  the  error  is  esti¬ 
mated  to  be  3^. 

Crosstalk  -  Since  it  is  necessary  to  generate  a  reference  signal 
corresponding  to  the  received  signal  frequency  at  the  receiving  site  there 
is  always  a  possibility  of  the  undesired  addition  or  subtraction  of  the 
reference  signal  from  the  received  signal.  This  may  occur  in  several 
ways  including  crosstalk  between  tape  recorder  heads,  poor  ground  returns 
particularly  in  the  front  end,  common  mode  coupling  in  power  supplies, 
and  magnetic  coupling  to  the  antenna.  For  our  purposes  any  or  all  of  the 


146 


above  will  be  considered  as  crosstalk  and  an  error  estimate  for  the  com¬ 
bined  effect  stated.  It  should  be  noted  that  a  well  designed  system  will 
have  low  crosstalk  and  therefore  it  will  be  difficult  to  measure  this  level. 
Measurements  of  crosstalk  on  the  systems  employed  in  these  tests  have 
only  been  possible  when  the  antenna  was  shorted  and  integration  times  of 
two  hours  were  employed. 

Crosstalk  level  has  been  observed  to  be  constant  for  a  particular 
piece  of  equipment,  but  not  necessarily  the  same  from  one  unit  to  another. 
Given  a  fixed  crosstalk  level  the  degree  of  degradation  will  be  greatest  at 
those  frequencies  and  sites  where  the  received  signal  is  lowest.  If  Hawaiian 
operation  at  45  Hz  is  chosen  as  an  example  of  worst  case  condition  the  sig¬ 
nal  estimate  may  be  in  error  by  as  much  as  1.  .  Since  the  received  signal 

phase  is  random  with  respect  to  the  reference  frequency  phase  this  repre¬ 
sents  the  worst  possible  or  peak  error,  and  will  vary  from  run  to  run. 
Contamination  From  Other  Transmitted  Fields  due  to  Antenna  Misalignment 

The  measurements  are  intended  to  measure  the  transverse  horizontal 
magnetic  field  component  H^  since  it  is  this  component  which  will  dominate 
as  the  range  from  the  transmitter  is  increased.  It  is  true,  however,  that 
both  a  radial  (Hq)  and  a  vertical  (H^)  component  are  present;  their  strength 
being  dependent  on  both  range  and  bearing  relative  to  the  transmitter  an¬ 
tenna  pattern.  If  the  receiving  antenna  is  aligned  perfectly  (aligned  along 
the  vertical,  and  pointing  precisely  at  the  transmitter)  then  the  H^  and  H^ 
components  will  not  influence  the  signal  estimates.  It  is  estimated  that 
an  angular  mis-alignment  of  ±  2°  relative  to  the  range  vector  between 
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transmitter  and  receiver  could  exist  at  given  site.  This  error  could  result 
from  various  causes,  among  them  unknown  local  magnetic  anomalies  which 
could  affect  compass  readings,  mis-reading  the  compass,  and  errors  in 
sighting  in  the  antenna.  An  error  in  the  vertical  alignment  of  the  loop 
would  most  likely  result  from  the  use  of  an  inaccurate  level,  or  a  lack  of 
care  in  antenna  erection.  This  error  is  believed  to  be  small  and  will  be 
ignored  henceforth. 

The  effect  of  small  antenna  misalignment  from  the  true  bearing 

will  introduce  an  error  into  the  estimate  of  H  since  the  H  and  H  fields 

cp  cp  p 

will  combine  vectorially.  Assuming  that  the  receiving  loop  antenna  is 
misaligned  by  0^  the  worst  case  resultant  voltage  induced  into  the  loop 
antenna  can  be  expressed  as  (see  equation  a-1) 

V  =  U)Nu  A  [  |h  I  cos  0  ±  I  H  I  sin  0  ]  (i-1) 

oc  o  1  cp  1  e  1  p  1  e 


where  H 

P 

direction. 


H  and  H 
cp  P 


is  the  magnetic  field  component  along  the  propagation  path 
For  the  single  horizontal  antenna,  the  far  field  ratio  between 
is  the  following  (Bannister  [18]  ) 


H 


-  1  =  1  — - —  tan  0 


(i-2) 


H 


cp 


UJdS 


where  0^  is  the  angle  between  electrical  axis  of  antenna  and  propaga¬ 

tion  direction 


is  the  complex  propagation  constant 
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ID  =  Z7t{ 


d  -  distance  from  receiver  to  transmitter 

The  errors  introduced  by  H  are  seen  to  greater  for  the  closer  in  site 

P 

locations.  Table  1-1  shows  what  we  believe  to  be  the  maximum  errors  in 
signal  estimation  due  to  antenna  misalignment. 

Table  1-1 


Maximum  Error  in  Signal  Estimation  due  to  Antenna  Misalignment 


— - 

75  Hz 

East- West 

East- West 

Antenna 

Antenna 

Utah 

1.  06* 

.  59* 

Hawaii 

.  27* 

.  15* 

It  is  concluded  that  the  maximum  error  at  any  site  presently  used  or 
contemplated  will  not  exceed  1.  06^  and  this  estimate  is  adopted  as  repre¬ 
sentative  of  this  class  of  error. 
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Discussion  of  Total  Error 


For  the  purpose  of  this  discussion  all  errors  are  expressed  as  per¬ 
centage  of  the  quality  in  question.  All  measurements  are  assumed  to  have 
Gaussian  distribution  about  their  mean  value.  Although  there  is  some 
justification  to  assume  that  the  stated  error  is  seldom  exceeded  and 
therefore  can  be  considered  to  represent  two  or  three  standard  deviations 
(2a,  3a)  we  will  employ  the  errors  as  representative  of  one  standard 
deviation.  It  should  be  pointed  out  that  those  types  of  error  which  are 
not  constant  from  measurement  to  measurement  will  be  reduced  by  inde¬ 
pendent  measurements.  For  these  errors  the  standard  deviation  will  be 
reduced  by  \Tn  where  N  is  the  number  of  independent  measurements. 

Based  on  a  single  measurement  at  a  single  frequency,  the  most 
probable  error  in  estimating  signal  strength  in  the  absence  of  noise  is 
±  5.  23$  (±  .  44  dB).  Additional  measurements  under  the  same  conditions 
will  act  to  reduce  the  error.  For  example,  if  each  of  the  reduceable 
sources  of  error  was  derived  after  3  independent  measurements  then  the 
probable  error  would  be  3.  21$  (±  .  27  dB).  In  actuality,  each  number 
derived  in  this  report  would  have  to  be  examined  as  to  the  number  of 
component  measurements  that  went  into  its  makeup,  an  error  estimated 
for  it,  and  then  combined  appropriately  with  other  estimates  with  which 
it  is  subsequently  used.  In  order  to  avoid  this  complexity  we  have  adopted 
the  course  of  assigning  a  worst  case  error  which  may  be  applied  to  all 
results. 
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Synopsis  of  Error  Budget 

Fixed  Sources  of  Error 

Turns  Area  .  8  $ 

Calibration  Pad  .  1  $ 

Antenna  Misalignment  1.  06$ 

Sub  Total  RMS  1.  33$ 

Reduceable  Sources  of  Error 

(Reduceable  by  V"N  where  N  =  no.  of  measurements) 

Calibration  Level  -  Field  1.0$ 

Frequency  Gain  Difference  2.  0  $ 

Calibration  Level  -  Laboratory  3.0$ 

Signal  Level  -  Laboratory  3.0$ 

Crosstalk  1.6$ 

Sub  Total  RMS  5.  06$ 

Total  RMS  Error  Based  on  One  Measurement  5.  23$ 

.  44  dB 

Probable  RMS  Error  After  N  =  6  2.  46$ 

.  21  dB 
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After  careful  consideration  of  the  total  number  of  measurements 


needed  to  obtain  an  estimate  at  a  single  frequency  it  was  concluded  that 
the  average  number  was  6  (N  =  6)  and  that  the  error  to  be  associated  with 
a  single  frequency  estimate  should  be  ±  2.  46$  (±  .  21  dB).  Of  course  all 
ancillary  estimates  which  make  use  of  more  than  one  single  frequency 
measurement  will  have  less  error  attached  to  them. 

An  example  of  the  use  of  the  error  estimate  as  it  pertains  to  the 
estimation  of  the  propagation  constant  follows.  Assuming  that  Utah  and 
Hawaii  are  the  two  sites  of  interest,  we  assume  that  the  signal  estimate 
in  Utah  is  high  by  .  Z  1  dB  while  the  corresponding  estimate  in  Hawaii  for 
that  measurement  at  that  frequency  is  low  by  .  21  dB.  This  leads  to  a 
differential  error  of  .  42  dB.  Noting  that  4.  52  megameters  separates  Utah 
and  Hawaii  we  conclude  that  the  error  in  estimating  the  propagation  con¬ 
stant  in  the  absence  of  noise  is 

—  4Z-  =  .  093  «  .  1  dB/Mm. 

4.  52 

Increasing  the  baseline  measurement  distance  will  reduce  the  error  due  to 
measurement  error  but  will  decrease  signal  strength  at  the  distant  station 
with  a  corresponding  increase  in  uncertainty  due  to  a  lower  S/N  -ratio. 
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